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1.  INTRODUCTION 


In  a  recent  article  Foutz  (1980)  introduced  a  new  test  for  goodness-of- 
fit,  to  be  called  the  test  in  the  sequel.  Although  the  test  was  proposed 
for  fitting  a  continuous  p-variate  distribution,  it  applies  equally  well  to 
univariate  problems.  The  null  distribution  of  the  test  statistic  was  shown 
to  be  distribution  free  as  well  as  being  independent  of  the  number  of  vari¬ 
ates  p.  Foutz  obtained  an  integral  representation  for  the  null  CDF  of  F  ; 

n 

explicit  expressions  for  this  CDF  were  given  for  sample  sizes  2  and  3. 

Closed  form  solutions  for  the  CDF  for  larger  sample  sizes  are  quite  hard 

to  derive  and  Foutz  has  provided  a  large  sample  normal  approximation  to  the 

null  distribution  of  F  . 

n 

In  a  preliminary  comparison  with  ten  replications  of  50  simulated  samples 
from  (I)  a  mixture  of  uniform  distributions  and  (Li)  a  standard  normal  distri¬ 
bution,  Foutz  found  that  the  F  test  outperformed  both  the  Chi  squared  test 

n 

and  the  Kolmogorov-Smirnov  (K-S)  test. 

In  this  paper  we  present  the  results  of  an  extensive  investigation  to 

compare  the  three  goodness-of-f it  tests,  the  Chi  squared  test,  the  K-S  test 

and  the  F  test.  Members  from  three  families  of  distributions,  viz.,  the 
n 

family  of  asymmetric  stable  dis tributions,  mixtures  of  normals,  and  the 
Pearson  family  have  been  selected  to  represent  the  true  underlying  distri¬ 
bution  of  the  samples.  T.he  goodness-of-f  it  tests  are  applied  to  test  the 
hypothesis  that  the  samples  are  from  a  standard  normal  distribution.  The 
measure  of  comparison  used  in  the  study  is  the  empirical  power,  based  on 
5000  replications  of  each  of  the  tests. 

The  Chi  squared  and  Kolmogorov-Smirnov  statistics  were  computed  using 
the  IMSL+  routines  GFIT  and  NKS1.  The  number  of  cells  used  in  the  Chi 

"^International  Mathematical  and  Statistical  Libraries,  7500  Bellaire 
Boulevard,  Houston,  TX  77036 
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squared  computation  was  6,  8,  and  8  for  sample  sizes  of  20,  30,  and  50, 
respectively. 

A  brief  discussion  of  the  F^  test  is  given  in  Section  2  and  a  descrip¬ 
tion  of  the  simulation  is  in  Section  3.  The  results  of  the  simulation  are 
presented  in  Section  4.  FORTRAN  codes  used  for  the  simulation  and  detailed 
tables  of  simulation  results  are  in  Appendices  I  and  II. 
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2.  F  TEST 
n 

The  F^  test  is  based  on  a  comparison  of  a  continuous  empirical  distri¬ 
bution  function  (CEDF)  with  the  hypothesized  CDF.  The  CEDF  is  obtained  by 
,fspreadingM  the  total  mass  over  "statistically  equivalent  blocks"  generated 
by  the  sample.  As  shown  in  Anderson  (1966)  and  Foutz  (1980),  given  the 
order  statistics  of  a  random  sample  of  size  (n-1) ,  n  statistically  equiva¬ 
lent  blocks  that  partition  the  sample  space  can  be  constructed  in  many 
different  ways  by  choosing  what  are  called  cutting  functions.  An  intuitively 
appealing  set  of  blocks,  which  is  the  one  used  in  this  study  is  obtained  by 
choosing  the  identity  functions  for  the  cutting  functions.  In  this  case, 
the  n  statistically  equivalent  blocks  are  =  (-<»,  x^],  ^  =  (x^^x^jl,.. 

B  =  (x,  n,°°)  where  x/ix,  , .  .  . ,  x.  are  the  order  statistics  of 

n  (n-1;  (1;  (2;  (n-1; 

a  sample  of  size  n-1.  The  CEDF  is  constructed  by  spreading  a  mass  —  cont- 

n 

tinuously  and  in  the  same  proportion  as  the  hypothesized  CDF  over  each 

block.  If  H  is  the  hypothesized  CDF  and  H  the  CEDF,  the  test  statistic 
o  n 

F  is  defined  as 
n 

F  =  Sup  |h  fx)  -  H  (x)  |  (1) 

n  xr  1  n  o  1 

Let  ,  i  =  1,2, ...  ,n  be  the  probability  contents  of  the  blocks  B_^ 

under  the  null  CDF  H,i.e.,D.  *  P[x  e  B.lH  ].  A  computationally  convenient 

O  1  1  o 

form  for  F  can  be  shown  to  be 
n 

n  1 

F  =  E  max(0, - D.)  .  (2) 

n  .  .  n  l 

i=l 

Foutz  has  provided  the  null  distribution  of  F  in  integral  form  and  derived 

n 

closed  form  solutions  for  n  =  3,4.  Simplifying  the  expressions  given  by 
Foutz  for  n  =  3,4  yields 


6x 

P(F3  <  x)  =  j  1  -  3(j  -  X)2 


0<x<  — 

J  <  x  <  2/3 
x  >  2/3 


(3) 


and 


20xJ 


P(F4  <  x) 


3  2  9  1 

-20x  +  18x  -  fx  +  tt 
4  16 

1  -  4(f  -  x)3 


0  <  x  <  7* 

-  -  4 

1  1 

—  C  X  c 

4  -  -  2 

*>!• 


(4) 


We  have  also  obtained  the  expression  for  the  null  CDF  for  n  =  5, 

4 


P(F5  <  x)  =< 


70x 

-105x4  +  80x2  -  12x2  +  -j|-x  -  Y25- 
45x4  .  80x3  +  2|8x2_176+_|I 

1  -  5(j  -  x)4 


0  <  x  <  — 

1  2 

T  <  x  -  T 

2  3 

5  X  -  5 

3  4 

5  <  X  -  5 

4 

X  >  5  * 


(5) 


As  is  evident  the  exact  distribution  is  quite  difficult  to  obtain  for 
higher  sample  sizes.  A  large  sample  normal  approximation,  due  to  Foutz  is 
given  by 


lim  P 
n-*” 


F  <  x 

n  - 


=  $ 


n(x  -  e 

-1  _2  l 

(2e  -  5e  ) 


(6) 


where  $  is  the  standard  normal  CDF.  For  n-1  =  20,30,50  we  used  this  approxi¬ 
mation  to  test  the  hypothesis  that  a  simulated  sample  from  U[0,1],  a  uniform 
distribution  on  [0,1],  is  in  fact  from  that  distribution.  In  80,000  repli¬ 
cations,  the  observed  significance  level  (number  of  hypothesis  rejections/ 
80,000)  was  consistently  smaller  than  the  nominal  value  as  can  be  seen  from 
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TABLE  1 


EMPIRICAL  SIGNIFICANCE  LEVEL  OF 
FOUTZ  Fn  TEST  USING  ASYMPTOTIC  APPROXIMATION 
(80,000  REPLICATIONS) 


Sample  Size 

20 

30 

50 

Significance  Level 

.10 

.0757 

.0800 

.0859 

.05 

.0372 

.0399 

.0428 

.01 

.0082 

.0083 

.0093 

9 


Table  1.  We  therefore  constructed  a  Monte  Carlo  CDF  of  F  for  n-1  =  2,3,4, 

n  9 

20,30,50  based  on  25,000  computer  generated  F  values;  these  represent  values 

n 

of  the  F^  statistic  for  testing  the  hypothesis  that  a  set  of  samples  from 

U[0,1]  is  in  fact  from  U [ 0 , 1 ] .  A  comparison  of  the  Monte  Carlo  CDF  with 

the  exact  CDF  for  n-1  =  2,3,4  is  provided  in  Table  2.  It  can  be  seen  that 

the  Monte  Carlo  CDF  provides  a  reasonable  approximation  even  for  small  n. 

The  power  properties  of  the  F^  test  detailed  in  this  paper  are  based  on 

the  Monte  Carlo  CDF  of  F  .  Critical  values  obtained  from  the  Monte  Carlo 

n 

simulation  for  significance  levels  .01,  .05,  .1  and  n-1  =  20,30,50  are  in 
Table  3.  In  Table  4  we  present  the  observed  significance  level  in  225,000 
replications  when  testing  if  a  set  of  samples  from  U[0,1]  is  in  fact  from 
U[0,1]. 
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TABLE  2 


MONTE  CARLO  SIMULATION  VS  EXACT  VALUES  OF  CDF 


n=3 

n=4 

n=5 

MC 

EXACT 

MC 

EXACT 

MC 

EXACT 

X 

.40 

.7883 

.7867 

.7624 

.7625 

.7607 

.7600 

.45 

.8604 

.8592 

.8725 

.8725 

.8671 

.8693 

.50 

.9182 

.9167 

.9396 

.9375 

.9405 

.9405 

.55 

.9598 

.9592 

.9694 

.9680 

.9768 

.9778 

.60 

.9864 

.9867 

.9860 

.9865 

.9915 

.9920 

.65 

.9989 

.9991 

.9952 

.9960 

.9975 

.9975 
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TABLE  3 


CRITICAL  VALUES  FOR  Fn  TEST 
OBTAINED  BY  MONTE  CARLO  SIMULATION 


Sample  Size 

20 

30 

50 

Significance  Level 

.10 

.42714 

.41903 

.40816 

.05 

.44865 

.43553 

.42116 

.01 

.48659 

.46579 

.44487 
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TABLE  4 


EMPIRICAL  SIGNIFICANCE  LEVEL  OF 
FOUTZ  Fn  TEST  USING  MONTE  CARLO  APPROXIMATION 
(225,000  REPLICATIONS) 


Sample  Size 

20 

30 

50 

Significance  Level 

.10 

.1006 

.0970 

.1003 

.05 

.0486 

.0486 

.0498 

.01 

.0103 

.0101 

.0102 
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3.  DESCRIPTION  OF  SIMULATION 


In  our  simulation  we  generated  deviates  from  three  families  of  distri¬ 
butions.  The  family  of  asymmetric  stable  distributions  has  been  previously 
used  by  Saniga  and  Miles  (1979)  to  investigate  the  power  of  several  goodness- 
of-fit  tests.  We  used  the  same  set  of  parameter  values,  a  =  1.0(. 3)1.9  and 
B  =  0(. 25)1.00  they  used.  Mixed  normal  distributions  have  often  been  used 
for  such  tests,  and  we  have  included  a  family  which  is  a  composite  of  N(0,1) 
and  N(0,a),  a  ^  1,  and  another  set  which  is  a  composite  of  N(.5,l)  and 
N(0,a).  Pearson  distributions  considered  include  a  variety  of  shapes,  from 
'hear”  normal  to  U  and  J  shaped.  Discussion  of  the  procedures  used  to 
generate  pseudorandom  deviates  from  each  of  these  families  follows. 

The  sample  data  was  obtained  by  starting  with  one  or  more  uniformly  dis¬ 
tributed  pseudorandom  deviates.  These  were  generated  using  the  IMSL  sub¬ 
routine  GGUBS,  the  basis  of  which  is  discussed  in  Lewis,  Goodman,  and 
Miller  (1969) . 

3.1  Generation  of  Asymmetric  Stable  Deviates  (Random  Stabilized  Standard 
Form) 

This  family  of  distributions  contains  as  a  special  case  the  normal 
distribution  (a  =  2 ,  B  =  0)  and  the  Cauchy  distribution  (a  =  1,  B  =  0) .  For 
a  <  2  the  distributions  have  infinite  variance,  which  makes  them  useful  for 
determining  the  ability  of  a  goodness-of-f it  test  to  detect  heavy  tailed 
distributions.  The  deviates  were  generated  using  the  program  given  in 
Chambers,  Mallows,  and  Stuck  (1976).  This  subroutine,  RSTAB,  used  one  deviate 
from  U [0 , 1 ]  and  one  exponentially  distributed  deviate  to  generate  one  RSSF 
deviate. 
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3.2  Generation  of  Mixed  Normal  Deviates 


Mixed  normals  of  the  form  (l-y)N(y^ ,a^)  +  y  NCy^jO^)  were  generated  using 

the  IMSL  routine  MDNRIS  to  convert  uniform  samples  into  standard  normal 

samples.  To  obtain  a  set  of  N  mixed  normal  variates  we  proceeded  as  follows: 

(i)  2N  uniform  random  variates  {u^}  were  generated  using  GGUBSr  (ii)  for 

each  i  =  1,...,N,  MDNRIS  was  used  to  convert  u.  to  a  standard  normal  z.. 

1  1 

was  then  transformed  to  a  normal  with  mean  y^  and  variance  a ^ ,  or  with 

mean  yn  and  variance  an,  depending  on  whether  u.fx  >  Y  ■>  or  u.1XT  <  Y, 

2  2  l+N  -  i+N 

respectively. 


3.3  Generation  of  Pearson  Type  I,  II  Deviates 

The  generation  of  samples  from  Pearson  Type  I  and  II  distributions  was 
done  via  table  look  up  and  linear  interpolation  on  the  inverse  cumulative 
distribution  function.  Sufficient  entries  to  assure  four  significant  decimal 
places  in  the  final  answer  was  achieved  adaptively  using  numerical  integration. 
Before  discussion  of  the  precise  details  of  the  process,  we  digress  for  a 
discussion  of  the  Pearson  distributions,  and  particularly  types  I  and  II. 

Following  Johnson  and  Kotz  (1970),  the  Pearson  probability  density  function 
p(x)  is  given  by 


1  d£ 
p  dx 


a  +  x 


2 

c  +  c,x  +  c0x 
o  1  2 


It  can  be  shown  that  a,  c  ,  X- ,  and  cn  can  be  expressed  in  terms  of  non- 

o  1  l 

negative  parameters  8  ,  and  t^ie  variance  y^,  obtaining 

Co  =  (4B2  ~  361)(10e2  -  1261  -  18)-1u2 

a  =  cx  =/f^  (82  +  3)(1082  -  12^  -  18)'1  /uj”  (7) 

c2  =  (282  -36-l  -6)(1082  -  12^  -  18)_1  . 
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Type  I  distributions  are  characterized  by 


12  -I 

K:  =  4  C1  (coC2)  <  0  ’ 

while  Type  II  have  k  =  0  with  6^  =  0,  <  3.  For  Type  I  and  II  the  prob¬ 

ability  density  function  is  of  the  form 


ml  m2 

p(x)  =  K(x  -  a1>  X(a2  -  x)  Z 


(8) 


Generation  of  a  sequence  of  pseudorandom  Pearson  Type  I  and  II  deviates 
involves  the  following  steps:  (i)  generate  a  sequence  of  deviates  {u_^}, 
from  U[0,1];  (ii)  transform  these  to  Pearson  deviates  by  finding  v.  so  that 

t Ti 

/  p(x)dx  =  u_^.  This  necessitates  being  able  to  efficiently  obtain  the 


■/ 


p(x)dx,  then  we  need  v_^  =  F  ^(u_^).  We 


inverse  CDF,  i.e.,  if  F(v)  = 

will  denote  F  ^  by  G  in  order  to  simplify  notation. 

Representation  of  the  inverse  CDF,  G(s),  was  achieved  by  linear  inter¬ 
polation  in  a  table  generated  by  numerical  integration  of  p(x).  The  adap¬ 
tive  process  used  to  assure  four  decimal  place  accuracy,  i.e.,  magnitude  of 

-4 

the  error  less  than  .5  x  10  is  described  now. 

The  error  in  linear  interpolation  between  points  (s^,G_^)  and  (s^_^,  ^i+l^ 

12  M 
is  no  more  than  —(s.+^  -  s^)  Q9  ,  where  G2  =  max  |g  (s)  |  . 

s . <s<s  _ 

1-  -  1+1 

Since  the  intervals  will  be  small,  we  can  approximate 

and  then  using  the  fact  that 


G2  by 


G,(si+i)  -G’(Si) 


S..T  -  S. 
1+1  1 


G'(s)  =  (F  1(s))’  =  rrr1 


F'(G)  P(G) 


,  we  obtain 
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P(1W  p<Gi> 

si+l  -  si 


Thus,  an  error  estimate  in  (s^,  is  given  by 

1  1 


p^Gi+l^  P(Gp 


s .  -  s  . 

1+1  1 


(SH-1  -  si>' 

8 


p(Gi+l)  p(Gi^ 


(s 


i+1 


Si} 


(9) 


Potential  problems  occur  if  p(G)  =  0 ,  as  may  happen  at  and  a^.  Since  we 

— 

require  error  less  than  .5  x  10  ,  it  is  clear  this  must  be  the  case  when 

i  i  —4 

-  Gj<.5  x  10  ;  hence  if  p(G)  is  very  small  in  the  interval 

we  have  an  alternative  scheme  for  accepting  an  interval,  one  which  came  into 
play  for  ranges  where  p  is  small. 

These  ideas  were  the  basis  for  adaptive  construction  of  a  suitable  table 


(s^,G^)  for  the  inverse  cumulative  distribution  function.  We  have 

s  =  0,  G  =  a.  .  We  describe  the  general  scheme  for  obtaining  (s .  ,  .  ,  G.,-) 
o  o  1  i+l  i+l 

given  and  note  how  an  estimate  of  G_^+^  is  generated  afterwards.  Given 

an  estimate  for  G.,-,  the  value  of  As.  =  s.(1  -  s.  is  obtained  by  numerical 
i+l  l  i+l  l 


numerical  integration  of 


/ 


Gi+i 


p(x)dx.  This  is  accomplished  by  the 


adaptive  quadrature  routine,  DCADRE,  from  the  IMSL  library.  If  m^  or  m^  are 

negative,  subtracting  out  the  singularity  was  used  for  intervals  near  a^  and 

a^,  respectively.  The  double  precision  version  of  DCADRE  is  used  and  an 

—6 

absolute  error  tolerance  of  10  is  requested.  The  routine  returns  As^  and 

-4 

the  error  estimate  (9)  is  computed.  If  it  is  less  than  .5  x  10  ,  the 


result  is  accepted,  we  set  s . , „  =  s.  +  As . ,  and  proceed  to  the  next  interval. 

i+l  l  i 

-4  2 

Suppose  the  error  estimate,  E  >  .  5  x  10  .  Since  E  =  0((AG)  ),  and 

est  est 
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O(AG^)  =  O(As^),  we  obtain  a  new  estimate  for  AG_^  by  taking  it  to  be 


.25  x  10 


est 


-4 


+  1 


This  takes  the  error  based  on  the  new  AG_^  to  approximately  midway  between  its 

—  /j. 

current  value  and  .25  x  10  .  More  than  one  correction  of  this  sort  may  be 

required,  in  particular  initially  where  a  reasonable  estimate  of  AG  is  not 

o 

available. 

Once  an  interval  has  been  accepted  (or  rather,  a  point  (s^Gj)  we 
increment  the  interval  counter  i,  and  then  estimate  the  new  interval  size 


AG.  by  AGi-l 
1  2 


.475  x  10 


est 


-4 


+  1  1.  This  yields,  based  on  the  above  assumptions, 


a  AG^  which  should  give  an  error  for  the  next  interval  which  is  midway  between 

—4 

the  previous  one  and  .475  x  10 

An  initial  value  for  AG  is  required.  Although  not  very  sophisticated, 

o 

a2  ”  ai 

we  simply  take  A Gq  =  —  >  and  depend  on  the  adaptive  machinery  described 

above  to  decrease  it  to  meet  the  error  tolerance,  or  increase  successive 
intervals  as  required  for  efficient  representation. 

It  is  true  that  the  final  value  of  s^,  call  it  s^,  should  be  equal  to  one. 
Because  of  the  numerical  integration,  this  is  never  exactly  the  case.  In 
the  worst  case,  3^  =  .01,  =  1-9,  we  obtained  the  final  value 


s._  -  .9999999845  an  error  of  about  1.55  x  10 
N 


-8 


In  order  to  avoid  any  prob¬ 
lems  due  to  the  table  not  covering  [0,1]  exactly  the  simplest  procedure 


was  to  replace  each  computed  s_^  by  s^/s  ,  thus  distributing  the  error  over 

the  entire  interval  and  yielding  a  consistent  table.  Note  that  this  is 

-4 

well  within  the  error  tolerance  of  .5  x  10 

The  procedure  has  been  thoroughly  tested  for  its  efficiency  in  repre¬ 
senting  the  inverse  CDF  as  well  as  for  accuracy.  For  the  particular  cases 
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in  which  we  were  interested,  the  inverse  CDF  was  represented  by  a  table  with 
no  more  than  729  entries.  This  occurred  for  the  case  6^  =  .25,  $2  ~  3.2, 
which  has  an  inverse  CDF  with  very  large  slopes.  For  more  gently  sloping 
inverse  CDFfs  we  were  able  to  use  as  few  as  101  intervals,  as  in  the  case 

=  .01,  £$2  “  1*75.  Most  intervals  had  an  error  estimate  of  between  .4  x  10 

-4 

and  .5  x  10  ,  which  shows  that  the  interval  sizing  process  we  have  used 

worked  quite  efficiently,  with  few  initial  estimates  being  rejected  for  being 
too  large,  without  also  resulting  in  intervals  much  too  small. 

The  routine  was  checked  against  the  published  tables  of  Johnson,  Nixon, 
Amos  and  Pearson  (1963)  for  many  values  of  the  parameters  and  at  most  of  the 
percentage  points.  With  a  few  exceptions,  where  a  difference  of  one  in  the 
fourth  decimal  place  was  noted,  the  results  check  exactly.  Generally  in 
these  cases  the  fifth  place  was  four  or  five  so  that  the  actual  error  was 
probably  well  within  our  tolerance. 


\ 
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4 .  RESULTS 


The  results  of  the  simulation  are  summarized  in  Tables  5-12.  The  empi¬ 
rical  power,  in  5000  replications,  of  the  Foutz  F  test  compared  with  that 

n 

for  the  Chi  squared  test  or  the  K-S  test  is  presented  as  a  percent  improve¬ 
ment  in  the  probability  of  rejecting  the  null  hypothesis  that  the  distri¬ 
bution  of  the  samples  is  the  standard  normal.  A  negative  entry  means  that 
the  power  of  the  F^  test  was  smaller  than  that  for  the  Chi  squared  test  or 
the  K-S  test,  whichever  is  appropriate. 

The  simulation  has  revealed  that  the  F  test  is  better  than  the  Chi 

n 

squared  test  which  in  turn  is  better  than  the  K-S  test  when  the  true  distri¬ 
bution  of  the  samples  is  heavy  tailed.  Many  such  distributions  are  included 
in  the  asymmetric  stable  family  as  well  as  the  family  of  mixtures  of  normals. 
For  the  mixed  normal  family,  if  the  two  normals  involved  in  the  mixture 
differ  in  the  means  the  K-S  test  performed  better  than  the  F^  test  even  when 
the  variances  differed.  We  now  discuss  in  more  detail  the  results  for  each 
of  the  three  families  of  distributions. 

4.1  Asymmetric  Stable  Family 

The  results  for  the  F  test  versus  the  Chi  squared  test  are  summarized 

n 

in  Table  5.  The  F  test  outperformed  the  Chi  squared  test  for  n  equal  to 
n 

21  and  31.  When  n  ■  51,  as  a  +  B  increased  the  performance  of  the  F^  test 
deteriorated  as  can  be  seen  from  the  lower  right  part  of  Table  5.  Another 
general  observation  is  that  as  the  significance  level  is  decreased  the  improve¬ 
ment  in  power  for  the  F^  test  is  accentuated. 

The  comparative  figures  for  the  F^  test  versus  the  K-S  test  are  pre¬ 
sented  in  Table  6.  Here  again  the  F  test  did  much  better  than  the  K-S  test; 
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the  results  also  indicate  that  for  the  asymmetric  stable  family  the  Chi 
squared  test  has  a  higher  power  than  the  K-S  test. 

4.2  Mixtures  of  Normals 

The  mixed  normal  distributions  that  we  considered  were  of  two  basic 

types.  The  first  type  is  of  the  form  (l-y)N(0,l)  +  yN(0,cr)  with  a  =  2,3,4 

and  y  =  .1,  .2,  .3,  1.0;  note  that  when  y  =  1.0  the  distribution  is  not  a 

true  mixture  but  N(0,a).  The  second  type  is  a  mixture  of  the  form 

(l-y)N(.5,l)  +  yN(0,a)  with  a  =  3,  y  =  .2,  .3  and  a  =  4,  y  =  .2. 

The  test  did  significantly  better  than  the  Chi  squared  test  except 

for  n  =  51  and  y  =  1.0  (see  Table  7);  in  the  latter  case  the  Chi  squared 

test  turned  out  to  be  the  better  of  the  two  tests. 

When  the  F  test  was  compared  with  the  K-S  test  the  F  test  turned  out 
n  n 

to  be  consistently  better  (Table  8) . 

In  the  case  of  a  mixture  of  the  second  type, which  included  a  location 

shift  (Tables  9,  10), the  K-S  test  proved  to  be  superior  to  the  F  test 

n 

while  the  comparison  between  the  F  test  and  the  Chi  squared  test  appeared 

n 

to  be  inconclusive. 

4.3  Pearson  Family 

We  chose  ten  distributions  of  types  I  and  II  to  encompass  a  variety 
of  shapes  as  shown  in  figure  1;  the  standard  normal  is  superimposed  as  a 
dotted  curve  in  each  of  the  graphs.  The  comparison  of  the  F^  test  versus 
the  Chi  squared  test  proved  to  be  inconclusive.  However,  the  K-S  test 
appeared  to  have  a  higher  power  than  the  F^  test  when  the  shape  of  the 
distribution  was  "near”  normal  such  as  the  ones  for 
(61,B2)  =  (0,2.3),  (0,2.8)  and  (.25,3.2). 
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TABLE  5 


n=21 


n-31 


n-51 


FOUTZ  F  TEST  VS  CHI  SQUARED  TEST 
n 

RANDOM  STABILIZED  STANDARD  DISTRIBUTIONS 


0-0 

0.25 

0.50 

0.75 

1.00 

Sig 

Level 

a 

1.0 

102.0 

102.7 

101.6 

95.8 

96.9 

.10 

1.3 

96.1 

88.6 

69.1 

49.1 

11.7 

1.6 

72.0 

68.3 

59.5 

49.0 

29.2 

1.9 

54.3 

55.4 

55.4 

53.9 

55.8 

1.0 

120.2 

117.2 

121.6 

114.9 

110.3 

.05 

1.3 

97.9 

98.2 

71.2 

49.6 

-  1.4 

1.6 

71.2 

63.4 

58.6 

41.4 

20.1 

1.9 

49.8 

45.6 

51.1 

49.8 

44.2 

1.0 

256.0 

215.9 

238.0 

234.5 

233.8 

.01 

1.3 

1.6 

174.8 

118.5 

158.6 

130.3 

131.8 

107.9 

71.7 

62.1 

-  8.0 

21.7 

1.9 

81.6 

73.0 

99.0 

89.7 

78.7 

1.0 

37.9 

36.8 

38.1 

39.3 

37.0 

.10 

1.3 

35.9 

34.0 

28.4 

13.8 

-10.4 

1.6 

26.6 

26.2 

16.5 

11.1 

0.9 

1.9 

17.6 

19.5 

19.4 

13.4 

16.7 

1.0 

47.4 

44.6 

46.6 

50.2 

46.7 

.05 

1.3 

41.8 

38.9 

31.5 

11.7 

-18.6 

1.6 

28.8 

25.7 

13.0 

5.2 

-  7.6 

1.9 

12.0 

15.8 

13.3 

8.4 

14.4 

1.0 

74.7 

75.2 

73.8 

73.5 

70.2 

.01 

1.3 

55.8 

46.2 

31.2 

3.9 

-33.7 

1.6 

32.2 

30.7 

11.0 

-  6.0 

-21.3 

1.9 

14.5 

18.1 

10.1 

8.0 

1.0 

1.0 

11.4 

13.3 

11.1 

10.6 

11.4 

.10 

1.3 

9.6 

12.3 

6.1 

-  1.1 

-11.1 

1.6 

1.2 

-  0.5 

-  4.9 

-  7.3 

-13.9 

1.9 

-10.8 

-  8.9- 

-11.0 

-11.1 

-11.2 

1.0 

17.4 

18.5 

17.1 

16.3 

17.1 

.05 

1.3 

12.6 

12.2 

5.0 

—3.5 

-18.2 

1.6 

0.5 

-  3.5 

-  8.7 

-12.0 

-21.7 

1.9 

-14.3 

-12.9 

-15.1 

-16.3 

-15.9 

1.0 

30.5 

30.5 

31.2 

28.4 

25.8 

.01 

1.3 

14.1 

13.8 

3.1 

-12.2 

-35.2 

1.6 

-  5.4 

-  8.3 

-20.6 

-23.5 

-37.0 

1.9 

-24.2 

-24.2 

-24.1 

-24.2 

-26.  L 
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TABLE  6 


FOUTZ  F  TEST  VS  KOLMOGOROV- SMIRNOV  TEST  RANDOM 
RANltf)M  STABILIZED  STANDARD  DISTRIBUTIONS 


Sie 

B-*  0.0 

0.25 

0.50 

0.75 

1.00 

Level. 

a 

1.0 

120.2 

123.8 

110.0 

120.4 

115.2 

* 

.10 

1.3 

102.8 

100.6 

76.0 

43.5 

0.9 

1.6 

87.7 

77.2 

61.3 

49.6 

22.7 

1.9 

65.5 

61.7 

62.9 

65.1 

57.8 

1.0 

197.5 

196.0 

186.7 

185.9 

192.1 

.05 

1.3 

170.8 

151.2 

111.3 

68.1 

-  2.8 

n=21 

1.6 

131.3 

110.4 

89.9 

67.0 

30.4 

1.9 

83.7 

90.5 

87.5 

85.6 

74.4 

1.0 

601.5 

465.5 

530.4 

522.7 

526.8 

.01 

1.3 

396.4 

338.1 

245.0 

121.4 

-  0.6 

1.6 

305.7 

243.9 

165.2 

117.7 

45.6 

1.9 

140.1 

147.5 

176.2 

176.5 

163.8 

1.0 

82.1 

80.6 

83.3 

84.7 

81.6 

.10 

1.3 

83.6 

77.6 

56.7 

25.5 

-  9.2 

1.6 

63.8 

59.5 

42.6 

29.4 

7.3 

1.9 

49.5 

55.6 

48.2 

44.8 

40.3 

1.0 

150.9 

149.8 

145.6 

152.7 

151.5 

.05 

1.3 

147.3 

129.2 

88.1 

40.5 

-  13.0 

n=31 

1.6 

107.8 

98.9 

65.0 

41.1 

11.4 

1.9 

77.2 

86.7 

69.5 

72.5 

68.2 

1.0 

428.5 

429.5 

447.0 

399.0 

471.3 

.01 

1.3 

386.1 

316.4 

190.6 

76.4 

-  16.2 

1.6 

267.3 

210.3 

136.9 

79.9 

24.1 

1.9 

190.8 

186.2 

153.0 

146.8 

136.8 

1.0 

38.5 

42.7 

38.9 

37.5 

37.7 

.10 

1.3 

49.3 

49.4 

29.6 

8.4 

-  10.3 

1.6 

40.0 

35.7 

24.0 

11.4 

-  6.7 

1.9 

24.4 

31.4 

26.1 

22.2 

19.5 

1.0 

77.6 

81.2 

78.0 

78.3 

76.0 

.05 

1.3 

94.4 

85.2 

53.2 

14.7 

-  16.4 

n=51 

1.6 

1.9 

79.4 

49.8 

61.4 

54.6 

37.9 

50.6 

17.0 

42.0 

-  8.1 

35.4 

1.0 

285.4 

287.0 

275.2 

270.9 

262.6 

.01 

1.3 

287.6 

243.0 

138.7 

37.7 

-  26.1 

1.6 

196.2 

173.3 

91.2 

41.0 

-  7.9 

1.9 

130.4 

126.3 

115.4 

111.8 

99.6 
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TABLE  7 


n=21 


n=31 


n=51 


FOUTZ  TEST  VS  CHI  SQUARED  TEST 
MIXED  NORMAL  DISTRIBUTIONS 


Y+  1.00 

0.30 

0.20 

0.10 

Sig 

N(0,a) 

Level 

- 

N(0,2) 

55.8 

69.3 

76.3 

40.7 

.10 

N(0,3) 

32.8 

86,8 

75.9 

44.7 

N(0 ,4) 

20.5 

101.4 

89.0 

75.7 

N(0,2) 

47.0 

46.4 

64.5 

23.7 

.05 

N(0 , 3) 

30.5 

71.3 

69.0 

34.7 

N(0,4) 

22.6 

91.1 

76.8 

48.1 

N(0,2) 

62.0 

262.1 

10  7.5 

88.6 

.01 

N(0,3) 

61.2 

109.8 

145.8 

86.3 

N(0,4) 

49.0 

118.3 

131.1 

125.0 

N(0,2) 

13.8 

46.9 

42.4 

22.0 

.10 

N(0, 3) 

4.3 

56.7 

50.6 

51.5 

N (0,4) 

1.8 

55.8 

58.7 

39.2 

N(0, 2) 

12.0 

38.3 

29.2 

28.6 

.05 

N(0 , 3) 

3.1 

53.2 

48.4 

53.2 

N(0,4) 

1.0 

57.4 

61.5 

44.4 

N(0 , 2) 

1.2 

62.1 

46.0 

52.1 

.01 

N(0 , 3) 

0.6 

47.3 

103.1 

104.7 

N(0,4) 

0.1 

84.0 

95.5 

61.5 

N(0,2) 

-13.1 

23.6 

26.1 

9.7 

.10 

N(0,3) 

-  5.9 

19.3 

28.4 

25.7 

N(0,4) 

-  1.4 

26.2 

29.5 

35.6 

N(0,2) 

-16.2 

28.9 

45.5 

17.4 

.05 

N(0 , 3) 

-  9.4 

29.2 

43.4 

38.0 

N(0, 4) 

-  2.5 

34.1 

45.4 

39.7 

N (0, 2) 

-26.9 

83.8 

54.8 

67.4 

.01 

N(0,3) 

-17.7 

41.6 

45.9 

67.9 

N(0,4) 

-  7.9 

24.5 

78.3 

51.4 
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TABLE  8 


FOUTZ  F  TEST  VS  KOLMOGOROV- SMIRNOV  TEST 
n 

MIXED  NORMAL  DISTRIBUTIONS 


T* 

1.00 

0.30 

0.20 

0.10 

_  Sig - 

_  N(0,a) 

Level 

N(0,2)  ' 

60.5 

36.1 

31.8 

8.7 

.10 

N(0 , 3) 

52.6 

59.3 

42.9 

19.3 

N(0 , 4) 

38.8 

71.5 

54.0 

40.2 

N(0 , 2) 

87.4 

49.5 

36.4 

6.8 

n=21 

.05 

N(0, 3) 
N(0,4) 

91.8 

74.4 

75.4 

90.2 

49.8 

66.6 

16.9 

43.7 

N(0,2) 

156.3 

94.4 

43.1 

15.8 

.01 

N(0,3) 

223.6 

112.3 

71.0 

41.8 

N(0,4) 

198.1 

107.0 

101.2 

26.6 

N(0,2) 

42.6 

19.3 

27.8 

13.9 

.10 

N(0,3) 

31.3 

52.9 

38.8 

24.1 

N(0,4) 

18.5 

66.2 

52.9 

21.1 

N(0,2) 

69.3 

29.9 

27.8 

20.6 

n=31 

.05 

N(0,3) 
N(0 ,4) 

60.8 

38.7 

65.7 

101.4 

52.9 

61.9 

36.8 

36.2 

N(0,2) 

135.8 

52.9 

67.3 

55.3 

.01 

N(0,3) 

183.9 

107.5 

120.3 

66.0 

N(0,4) 

129.8 

163.9 

152.2 

72.1 

N(0,2) 

22.4 

26.8 

23.6 

9.7 

.10 

N(0, 3) 

9.6 

47.8 

39.4 

21.2 

N(0,4) 

3.1 

66.4 

51.1 

36.3 

N(0,2) 

45.7 

30.0 

39.9 

9.0 

n=51 

.05 

N(0, 3) 
N(0,4) 

24.5 

10.6 

74.0 

110.7 

52.2 

73.3 

29.9 

52.3 

N(0 , 2) 

117.6 

61.9 

81.1 

20.3 

.01 

N(0,3) 

96.0 

215.7 

83.3 

39.1 

N(0,4) 

53.0 

l 

170.9 

164.5 

86.7 
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TABLE  9 


FOUTZ  F^  TEST  VS  CHI  SQUARED  TEST 
MIXED  NORMAL  DISTRIBUTIONS 


SIGNIFICANCE  LEVEL 


MIXED  NORMAL 

.10 

.05 

.01 

; 

0 . 7GxN (0.5, l)+0 . 30xN(0 .0,3) 

24.0 

-  0.7 

-20.8 

n=21 

0 . 80xN(0 . 5 , l)+0 . 20xN(0 .0 , 3) 

14.4 

-  9.2 

-27.2 

0 . 80xN (0 . 5 , 1) +0 . 20*N (0.0,4) 

17.2 

-  3.2 

-21.1 

0.70xN(0.5,l)+0.30XN(0.0,3) 

9.6 

-  6.8 

-28.6 

3 

II 

u> 

h-* 

0 . 80xN (0 . 5 , 1) +0 . 20xN (0.0,3) 

3.3 

-12.4 

-33.3 

0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,4) 

3.7 

-10.7 

-30.2 

0 . 70xN (0 . 5 , l)+0 . 30xN (0.0,3) 

25.1 

-24.1 

-43.6 

n=51 

0 . 80xN (0 . 5 , l)+0 . 2GxN(0 .0,3) 

-  2.3 

-23.6 

-48.6 

0. 80xN(0 . 5 ,l)+0 . 20xN(0. 0,4) 

-  4.1 

-26.3 

-45.5 

i 
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TABLE  10 


FOUTZ  F  TEST  VS  KOLMOGOROV- SMIRNOV  TEST 
n 

MIXED  NORMAL  DISTRIBUTIONS 


SIGNIFICANCE  LEVEL 


MIXED  NORMAL 

.10 

.05 

.01 

0. 70xN(0. 5 ,l)+0. 30xN (0.0,3) 

-17.0 

-28.9 

-33.7 

n=21 

0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,3) 

-27.3 

-38.2 

-40.7 

0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,4) 

-25.5 

-32.5 

-35.1 

0. 70xN(0. 5 , l)+0. 30xN (0.0,3) 

-25.4 

-34.5 

-43.4 

n=31 

0. 80xN (0.5, l)+0. 20xN(0. 0, 3) 

-34.8 

-45.5 

-50.7 

0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,4) 

-34.6 

-40.2 

-45.9 

0. 70xN(0.5,l)+0.30xN(0.0,3) 

-28.7 

-47.6 

-56.7 

n=51 

0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,3) 

-48.7 

-57.4 

-64.9 

0 . 80XN (0 . 5 , l)+0 . 20xN (0.0,4) 

-43.3 

-51.3 

-61.0 
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TABLE  11 


FOUTZ  F  TEST  VS  CHI 
n 

SQUARED  TEST 

PEARSON  DISTRIBUTIONS 

SIGNIFICANCE 

LEVEL 

'l 

S2 

.10 

.05 

.01 

0.0 

2.30 

3.4 

-10.5 

-20.6 

0.0 

2.80 

21.1 

-  0.9 

-12.4 

0.01 

1.75 

12.8 

41.1 

57.8 

0.01 

1.90 

-  9.1 

-  2.4 

18.9 

n=21 

0.25 

3.20 

16.6 

-  0.4 

-  2.3 

0.50 

3.00 

17.5 

5.8 

-  9.3 

1.00 

3.40 

-12.5 

-16.2 

-15.8 

1.00 

3.60 

9.1 

-31.4 

-30.9 

1.00 

3.80 

16.1 

-17.6 

-24.8 

2.00 

5.50 

-16.1 

-14.9 

-20.7 

0.0 

2.30 

-  6.6 

-19.0 

-22.5 

0.0 

2.80 

10.7 

-12.1 

-  5.4 

0.01 

1.75 

-  5.3 

39.4 

83.5 

0.01 

1.90 

-18.7 

-  6.9 

26.9 

n=31 

0.25 

3.20 

0.7 

-11.8 

-  3.6 

0.50 

3.00 

-  1.1 

-  0.6 

-12.2 

1.00 

3.40 

-25.5 

-23.5 

-23.4 

1.00 

3.60 

-  5.7 

-41.4 

-41.7 

1.00 

3.80 

0.0 

-26.9 

-32.9 

2.00 

5.50 

-29.5 

-19.4 

-26.5 

0.0 

2.30 

30.4 

-39.0 

-19.7 

0.0 

2.80 

34.9 

-40.0 

-17.3 

0.01 

1.75 

6.6 

68.1 

132.9 

0.01 

1.90 

-11.8 

-15.8 

20.7 

n=51 

0.25 

3.20 

33.3 

-10.8 

-15.6 

0.50 

3.00 

2.5 

-19.9 

-20.2 

1.00 

3.40 

-37.4 

-35.7 

-35.6 

1.00 

3.60 

-13.6 

-53.7 

-56.7 

1.00 

3.80 

-  4.7 

-40.6 

-46.5 

2.00 

5.50 

-38.5 

-28.7 

-38.0 
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TABLE  12 


FOUTZ  F  TEST  VS  KOLMOGOROV- SMIRNOV  TEST 
a 

PEARSON  DISTRIBUTIONS 


SIGNIFICANCE  LEVEL 


*1 

*2 

.10 

.05 

.01 

0.0 

2.30 

-22.8 

-25.0 

-26.2 

0.0 

2.80 

5.4 

-18.8 

-  7.2 

0.01 

1.75 

-  9.4 

20.6 

66.0 

0.01 

1.90 

-33.3 

-15.3 

16.5 

n=21 

0.25 

3.20 

-  8.1 

-14.8 

-  4.1 

0.50 

3.00 

9.9 

21.7 

34.9 

1.00 

3.40 

84.8 

94.4 

95.9 

1.00 

3.60 

52.9 

64.5 

76.7 

1.00 

3.80 

38.4 

44.0 

59.4 

2.00 

5.50 

62.4 

74.7 

70.4 

0.0 

2.30 

-28.9 

-36.8 

-36.9 

0.0 

2.80 

11.2 

-20.8 

-13.0 

0.01 

1.75 

-14.1 

27.6 

89.1 

0.01 

1.90 

-36.7 

-15.4 

19.4 

n=31 

0.25 

3.20 

-12.9 

-14.0 

-  9.3 

0.50 

3.00 

-  0.6 

25.0 

31.3 

1.00 

3.40 

100.0 

134.3 

131.7 

1.00 

3.60 

55.5 

79.4 

91.8 

1.00 

3.80 

39.9 

46.2 

65.3 

2.00 

5.50 

63.3 

91.0 

100.0 

0.0 

2.30 

-  9.1 

-48.6 

-40.0 

0.0 

2.80 

-  4.9 

-38.9 

-35.8 

0.01 

1.75 

-18.4 

62.0 

171.8 

0.01 

1.90 

-37.4 

-22.9 

26.6 

n-51 

0.25 

3.20 

-  8.1 

-  9.6 

-13.8 

0.50 

3.00 

0.0 

25.9 

47.0 

1.00 

3.40 

169.9 

225.8 

245.5 

1.00 

3.60 

64.9 

138.9 

154.4 

1.00 

3.80 

42.1 

69.0 

85.2 

2.00 

5.50 

103.0 

151.0 

164.7 

30 


5.  CONCLUDING  REMARKS 


The  superior  performance  of  the  Foutz  test  in  detecting  certain  types 
of  deviations  from  the  hypothesized  distribution  leads  to  several  more  prob¬ 
lems  to  be  considered.  Of  primary  importance  is  the  generation  of  percentage 
points  for  the  distribution  of  F^  for  various  values  of  n.  The  intractability 
of  the  problem  of  obtaining  the  exact  distribution  requires  an  empirical 
approach  to  finding  a  correction  to  the  asymptotic  approximation  given  by 
Foutz . 

Since  the  test  is  also  applicable  to  p-variate  distributions,  an  investi¬ 
gation  of  ways  to  obtain  the  statistically  equivalent  blocks,  and  then  the 
probability  content  of  them,  at  least  for  p  =  2,  is  to  be  considered.  The 
problems  of  obtaining  these  blocks  and  their  contents  becomes  increasingly 
complicated  in  higher  dimensions. 

The  authors  wish  to  acknowledge  extensive  use  of  the  IBM  360/67  computer 
system  at  the  Naval  Postgraduate  School,  on  which  all  computations  were 
performed. 
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APPENDIX  1 


this  PROGRAM  GENERATES  DEVIATES  tn  random  stabilized  standard  form 
BLGCKS  GF  20,  30  »  AND  50  DEVIATES  ARE  CONSIDERED  BY  GENERATING  50, 
THEN  CONSIDERING  THE  FIRST  20,  THE  NEXT  30,  AMD  ALL  50. 

THE  HYPOTHESIS  THAT  THEY  CAME  FROM  A  NORMAL  DISTRIBUTION  WITH 
MEAN  AND  VARIANCE  AS  GIVEN  IN  INPUT  IS  TESTED. 

THE  TESTS  PERFORMEC  ARE  THE  CHI  SQUARED,  KOLMOGOROV- SMI RNOV ,  AND 
FOUTZ  TEST  AT  THE  CCNFIOENCE  LEVELS  OF  10,  5,  AND  1  PERCENT. 

THE  INPUT  VALUES  APE  MEAN,  VAR! ANCE, I  SEED ,  AND  THE  NUMBER  OF  REPLI 
CATIONS. 

IMPLICIT  R£AL*8  (A-Z) 

INTEGER**  T,  IER»N»I»NP1,N1»N2»N3*I SEED, NR , N STOP , IDF  »IQ,ISEED1 

1  , NCEL ( 3 )  , NST ( 3 ) ,N SM P (3 i ,NCHI (3 ,3 ) ,NCHIN (3,3 ) , NKS (3, 3  )  , 

2  NKSN(2»3)»NFNN(2,3),NFN(3»3) ,NSR 
DIMENSION  RD(55  J  ,CHISQU(2)  , 

1  KOLSMR ( 2 ) ,C0UTZL<2)  ,FNTST(3,2)  ,  ZA LF ( 3 ) , RN (55  I 
REAL**  CELLS (50 l » CCM P ( 50 ) »  Q» CH I , R ( 55) , PD I F ( 6) , QTST( 3 ) , RNN( 55) , 

1  STST( 3  »  3 ) 

DATA  STST/. 26*73,. 21 756,. 165 5S,. 29*08, .2*170,. 188*1 ,.352*1, 

1  .28987,. 2260*/ 

DATA  NCEL/8,10,10/,NST/l,21,l/,NSMP/20,30,50/,QTST/.l,  .05,  .01/ 

1  , CHIS  QU / • CH I  SQUA ' , *  RED  TE ST • /,KOLSMR/ • KOLMOG-S • , • MI R  TEST*/ 

2  ,FOUTZL/' F0UT2  T E ST •/, ZALC1 , ZALC2/ . 2*306900, . 3678 79DQ / 

DATA  WT1,MEAN1 , VAR1 , WT2 ,ME AN 2 , VAR2 /. 8D0 ,0. DO , 1 .DO , .2  DO ,0 . DO , 3 . DO/ 

3  , ZAL F/ 1. 281 55C0, 1.6**8500, 2.3263500/ 


THE  NEXT  *  CARDS  PUT 
FOUTZ  DISTRIBUTION. 


IN  THE  EMPIRICAL  CRITICAL  VALUES  FOR  THE 


1 

2 


DIMENSION  FN=mp<3,3) 

DATA  FNEMP /. *271*00,. *190300,. *)816D0, 
.**865 CO , .*355300, .*211600,  .*865900, 

•  *6579 DO,.  ***870  0/ 


100 


1* 

1 


17 

12 


110 

5 

3 


REAL**  ALPHA, BETA ,BFRIME 
COMMON / NCR PRM/ ME AN  » VARI »  SVAR I 
EXTERNAL  NORM , UN  1 F 
T  =  1 

PIB2  =  CAT  AN(1.DC )*2.00 
PRINT  1 

READ  2, MEAN, VARI 
IF (MEAN+VARI .LE.O.CO)STOP 
FORMAT ( 2E5  .0 ) 

SVAR  I  *  DSORT  (VARI) 

WRITE ( * , 1* ) 

FORMAT ( 1H1 ) 

FORMAT  ( '+INPUT  MEAN  AND  VARIANCE  FORMAT(  2E5. 0  )  *  ) 

PRINT  17 

READ  12, ALPHA, BETA 

FORMAT  (  '  ♦  INPUT  ALPHA  ANO  BETA,  FORMAT!  2F5.0)  '  ) 

FORMAT (  5F5 .0  ) 

PP  INT  3 

READ  5 , ISEED»NB 

ISEED1  =  I  SEED 

BPRIME  =  BETA 

IF ( ALPHA .EC.O. )GC  TC  110 

PHIZ  =  -PI B2*BETA*(1.  -  ABS ( 1 •  -  ALPHA )) /ALPHA 
BPRIME  =  DTAN(PIB2*(1.  -  ALPHA  )  )*0TAN( PHI Z*ALPHA) 

CONTINUE 

DSEED  =  DFLOAT ( I  SEED ) 

FORMAT  (110,15) 

FORMAT ( • +INPUT  SEED  AND  NUMBER  OF  REPLICATIONS,  FORMAT! 110,15)* ) 

CO  205  IQ= 1 » 3 

DO  205  Nl*l,3 

NCHI(N1,IQ)  =  C 

NCHIN(NltlQ)  =  0 

NKS(N1 , I C)  =  0 

NKSN ( N1 , IQ )  =  C 
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0 


c 

c 

c 

c 


THE  FOLLOWING  CARD  PUTS  IN  THE  EMPIRICAL  VALUES  FOR  THE 


FOUTZ  TEST 


FNTST  (Nl , IQ )  =  FNEMP ( Nl, IQ ) 


C 


205  NFNN ( N 1 , IQ  )  =  0 
DO  300  1=1, NR 

CALL  RA NSRV( 50 ,R N  ,RD , CSEED , ALPHA , 8 PRIME ) 


N3  =  NSMP(Nl) 


CALL  GF IT (UN IF ,NC EL (Nl ) ,PNN , N3 , CELLS , COMP , CHI  , IDF,Q, IER) 

DO  241  I C= 1 , 3 

I F  ( Q  •  L  T •  GT  ST  ( 1 0  ) )NCHI(N1»IQ)  =NCHI(N1,IQ)  +  1 

241  CONTINUE 

CALL  NKS1(  UN  IF  ,RNN  ,  N3  ,  PDI F  ,  IS  R  ) 

CO  242  IC  =  1,3 

IF(PQIF(1)  .GT.STST(M,IQ))NKS(N1,  IC)  =  NKS(N1,IQ)  +  1 

242  CONTINUE 

CALL  FGUTZ(UNIF,PNN»N3,Q) 

DO  243  I C=1 , 3 

IF(0.GT.FNTST(N1, 10) )NFN(N1, IQ)  =  NFN(N1,I0)  +  1 

243  CONTINUE 
IDF  =  0 

CALL  GFIT(N0RM,NCEL(N1),R,N3, CELLS, COMP, CHI ,IDF,Q,IEP ) 

DO  244  IC  =  1,3 

IF(0«LT.QTST(IQI )NCHIN(N1,IQ)  =  NCHIN(N1,I0)  ♦  1 

244  CONTINUE 

CALL  NKS1(N0RM,R,N3,P0!F,I5R) 

DO  245  10=1,3 

IF(PDIF(1).GT.STST(M,IQ))NKSN(N1,IQ)  =  NKSN(N1,I0)  ♦  1 

245  CONTINUE 

CALL  FOUTZ (NCR N,R,N2»Q) 

DO  246  I C=1 , 3 

IF(O.GT.FNTST( Nl ,IC) )NFNN(N1 , IC)  =  NFNN(N1,I0)  +  1 

246  CONTINUE 
250  CONTINUE 
300  CONTINUE 

ISEED  =  CSEED 

WRITS (4, 7) MEAN,VAFI,NR,!SFED1, ISEEC 
WRITE! 4, 11 )ALP HA, BETA 

11  FORMAT ( // '  RANDOM  STABLE  STANDARDIZED  FORM  (RSSF),  ALPHA  =  ',F5.2, 
1  *,  BETA  5',F5.2//I 
WRITE  (4,  6KHIS0U,  ( ( NCHI (Nl , 1 0 )  , 

1  Nl  =  l,3)  ,10=1,3)  ,  ( ( NCHIN ( Nl , IC), 

2  Nl=l, 3)  ,10=1,3) 

WRITE(4,8)K0LSMR,( (NKS(N1, IQ),N1=1,3) ,10=1,3) , 

1  (  <NK$N(  M,IC)  ,N1  =  1,3)  ,10=1 ,3  ) 

WRITE! 4,8) FOUTZL, ( (NFN(N1, 10) ,N1  =  1  ,3) ,10=1,3 ) , (( NFNN (Nl, IQ) ,N1  =  1, 
1  ),  10=1,3) 

7  FORMAT  (A/*  THE  VALLES  OF  MEAN  (U)  AND  VARIANCE  (V)  ARE  S2F12.5/// 
9  *  THE  RESULTS  OF  ',15,'  P5PL ICATI CNS,  STARTING  WITH  SEED', 

1  I12//33X,  'ENDING  WITH  NEXT  SEEC',112//) 

8  FORMAT( ///20X,2A8//20X, '20  PIS', 

2  6X ,  '  30  PTS '  ,6  X ,  '50  PTS  V/10X,  '  10*  ',  31  12/'  CONTROL  ',  2X  ,  •  5?*, 
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END 


200 


SUBROUTINE  RANSR V ( NO  , RAN, R , D SEED , ALPHA » BPRIME ) 
IMPLICIT  REAL*8  (A-M, 0-Z) 

DI  MENSICN  R ( 1 )  ,RAN  <11 

REAL*4  RN<  100 >  .ALPHA .  BPRIME  ,R  STAB  ,  Vs 

CALL  GGUBS(DSEECt2*ND.RN) 

DO  200  NN= 1 »  ND 

PN<NN+NDI  ) 

R(NN)  =  RSTAB< ALPHA , BPRIME ,RN(NN) ,W> 

CONTINUE 

RETURN 

END 


SUBROUTINE  NCRN(X.P) 

REAL*8  MEAN.VAPI .SVARI 
COMMON /NCR  PR M/ ME  AN tVARI .SVARI 
T  =  (X  -  MEAN  I /S VAR  I 
P  =  .5*£RFC(-T*.7G71068) 
RETURN 
END 


subroutine 

P  =  X 

RETURN 

END 


UNIF(X,P) 


FUNCTION  RSTA6(ALPHA»8PRIMEtUtW) 

RSTAB  IS  A  RANDOM  STABLE  STANCARCI ZED  FORM 


(WHATEVER ) 


ARGUMENTS 
ALPHA  = 
BPRIME  = 
U  = 

W 


50 


CHARACTERISTIC  £XPCN£NT 
SKEWNESS  IN  REVISED  PARAMETERIZATION 
UNIFORM  VARIATE  ON  (0*1) 

EXPONENTIALLY  DISTRIBUTED  VARIATE 
DOUBLE  PRECISION  CA,CB 

DATA  PI BY2/1. 57079633/, PI8Y4/ .785398 163/ *THR1/. 99/ 

EPS  =  1.  -  ALPHA 
COMPUTE  SOME  TANGENTS 
PHIBY2  =  PIBY2*<U  -  .5) 

A  =  PHIBV2*TAN2(PHI8Y2) 

B8  =  TAN2(EPS*PH!EY2 ) 

B  =  EP  S*  PHI BY2*BE 

IF(EPS .GT.-.99ITAU  =  BPRIME /( TAN2 ( EP S*PI BY2 ) *P IBY2 ) 
IF(EPS.LE.-.99)TAL  =  BPRIME*P  I  BY2*5PS*(  1.  -E°S ) *TAN2< < 1. 
1  PIBY2 ) 

COMPUTE  SOME  NECESSARY  SUBEXPRESSIONS 
IF  PHI  NEAR  PI  BY  2*  USE  DOUBLE  PRECISION. 

IF(A.GT.THR1 )G0  TO  50 
SINGLE  PRECISION 
A2  =  A**2 
A2P  =  1.  +  A2 
A2  =  1.  -  A2 

B2  =  8**2 

B2P  =  1.  +  B2 
B2  =1.  -  82 
GO  TO  IOC 
DOUBLE  PRECISION 
DA  =  CBLE(A)**2 
DB  =  DBLE( B )  **2 
A2  =  l.CO  -  DA 
A2P  =  1  •  DO  +  DA 
B2  »  1.  -  OB 
D2P  =  1.  +  DB 
COMPUTE  COEFFICIENT 


-  EPS)* 
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100  Z  =  A2P* ( B  2  ♦  2.*PHIBY2*8B*TAU1/(W*A2*B2P) 

C  COMPUTE  THE  EXPONENTIAL-TYPE  EXPRESSION 

ALOGZ  =  ALOG(Z) 

0  =  D2(EPS*AL0GZ/(1.  -  EPS ) )*( ALOGZ/( 1.  -  EPS)) 

C  COMPUTE  STABLE 

P.STAB  =  (1.  +  EP$*D)*2.*((A  -  e)*(l.  +  A*B)  -  PHI  6Y2  *T  AU*BB*(  B*A2 
1  -  2  .*A  )  )/ (A2*E2P )  +  TAU*0 
RETURN 
END 


100 


FUNCTION  D2(  Z ) 

EVALUATF  ( EXP ( X  I  -  1)/X 

DOUBLE  PRECISION  F1,P2»G1,Q2,C3,PV,ZZ 

DATA  P1»P2»Q1,Q2»Q3/.  84006685 2 5364 €3 23 903 ♦•23001114158996456902 ♦ 
.  . 168013370507926648 C4, .1800 133704C7390023D3, 1 .DO/ 

THE  APPROXIMATION  1801  FOP  HART  ET  AL  (1968, P213) 

IF ( ABS (Z).GT .0.1  ICC  TO  100 

ZZ  =  Z*Z 

PV  =  PI  4  ZZ*P  2 

C2  =  2 . CO*PV/ ( 01  4  Z Z*(Q2  +  ZZ*Q3)  -  Z*PV) 

RETURN 

D2  =  (EXP(Z)  -  l.)/Z 

RETURN 

END 


FUNCTION  TAN(XARG) 

C  TANGENT  FUNCTION 

LOGICAL  NEG.INV 

DATA  P0,PltP2,00,Ql ,02/  . 129  22103553 ,-.887662 3 77 El, .528644456E-1 , 
1 .1645 2 93 32 £3,- .451*2056  IE  2, 1./ 

C  THE  APPRCX IMATI CN  4283  FROM  HART  ET  AL  (1968,  P.  251) 

DATA  PIBY4  / . 785 39 8 1 6 3/ , PI B Y 2 / 1 . 57C79633 / , PI /3 . 1 4159265/ 

NEG  =  .FALSE. 

INV  =  .FALSE . 

X  =  XARG 
NEG  =  X.LT.Q. 

X  =  ABS(  X) 

C  PERFORM  RANGE  RECUCTION  IF  NECESSARY 

IF(X.LS.PIBY4)GC  TC  50 
X  =  AMOD(X.PI) 

IF(X.LE.FIBY2)GQ  TO  30 
NEG  =  .NCT.NEG 
X  =  PI  -  X 

30  IF(X.LE. FI  BY 4)  GC  TC  50 

INV  =  .TRUE. 

X  =  PIBY2  -  X 
50  X  =  X/PieY4 

C  CONVERT  TO  RANGE  CF  RATIONAL 

XX  =  x  *x 

TAN  =  X*(PO  4  XX*(F1  4  XX*P2I)/(C0  4  XX*<01  +  XX*Q2») 

IF(NEG)TAN  =  -TAN 
IF(INV)TAN  =  1  ./TAN 
RETURN 
END 


FUNCTION  T AN2 (  XAR G  ) 

C  COMPUTE  TAN(X)/X 

C  FUNCTION  DEFINED  CNLY  FCR  AB S ( XARG ) . LE  .  P I  BY  4 

C  FOR  OTHER  ARGUMENT  RETURNS  TAN(X)/X,  COMPUTED  OIRECTLY 

DATA  P0,P1  ,P2, CO, Cl, C2/. 129221 03 5E3, -.887662 377Elt .528644456E- 1 , 
1  .16452523 2c 3, -.45122056 15 2,1. / 

C  THE  APPRCX IMATICN  4283  FROM  HART  ET  AL  (1968,  P.  251) 

DATA  PIBY4,PIEY2,PI / .7853981 6 3 , 1 . 57079633 , 3 . 14 15926 5 / 

X  «  ABS ( XARG ) 

IF(X.GT.PIBY4)GC  TC  200 
X  =  X/PIBY4 

C  CONVERT  TO  RANGE  CF  RATIONAL 
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THIS  PRCGRAM  GENERATES  MIXED  NORMAL  DEVIATES  WT1*N ( MEAN  1 , VAR  1 )  + 
( 1  -  WT1)*N!M£AN2  ,VAR2)  . 

BLOCKS  CF  20,  30,  AND  50  DFEVIATES  ARE  CONSIDERED  BY  GENERATING 
50,  THEN  CONSIDERING  THE  FIRST  20,  THE  NEXT  30,  AND  ALL  50. 

THE  HYPOTHESIS  THAT  THEY  CAME  FROM  A  NORMAL  DISTRIBUTION  WITH 
MEAN  ANO  VARIANCE  AS  GIVEN  IN  INPUT  IS  TESTED. 

THE  TESTS  PE RF CP. NEC  ARE  THE  CHI  SQUARED,  KOLMOGOROV-SMIRNOV,  AND 
FOUTZ  TEST  AT  THE  CCNFIOFNCF  LEVELS  OF  10,  5,  AND  I  PERCENT. 

THE  INPUT  VALUES  APE  MEAN,  VAR  IANCE, ISEED,  AND  THE  NUMBER  OF  REPLI 
CATIONS. 

IMPLICIT  REAL* 8  !A-Z) 

INTEGER**  T.IER.N.I  ,NP1 ,N1 ,N2  ,N3,I SEEO.NR.NSTOP, IDF, IC, ISEED1 

1  , NCEL ( 3  ) » NST !  3  ) »NSMP!3i*NCHI (3,3) ,NCHIN ! 3, 3) , NK S! 3 , 3 ) , 

2  NKSN<  3,3) ,NFNN( 3,3) »NFN (3 » 3  )  , NSR 
DIMENSION  RD(55) ,CHISQU(2) , 

1  K0LSMR(2) ,FQUTZL(2) , FNTST ! 3 , 3 ) , ZALF ! 3 ) , RN t 55 ) 

REAL*4  CELLS<  50  )  ,CCM F < 50) , Q, CH I , R ! 55  ) , PO I F ! 6 ) , 

I  STST ( 3 , 3 ) 

DATA  STST/. 26473, .21756, .16959, .29*08, .2*170, 

1  . 28987, .226C*/ 

DATA  NCEL/ 8, 10 , 10/ ,N ST/1, 2 1 , 1 /  ,N$MP/20 , 30 , 50/ , QTST/. 1 ,.05,. 01/ 

1  , CHI  SQL/' CHI  SCUA* ,'RE0  TEST  • /, KOLSMR/ • KOLMOG-S • , *MIR  TEST'/ 

2  » FOUTZ  L/*  FOUTZ  TE ST •/, ZALC1 , ZALC2/ .243069D0, .36787900/ 

DATA  WT 1, MEAN1, V AR 1, WT2, MEAN 2, VAP 2/. 800,0.00,1. DO,. 200,0.00, 3. DO/ 

3  ,  ZALF/ 1. 281 55 CO,  1.6**8  500, 2. 3  263 5D0/ 


QTST ( 3  ) , RNN (55) 
188*1, .352*1, 


THE  NEXT  *  CARDS  PUT  IN 
FCUTZ  TEST. 


THE  EMPIRICAL  CRITICAL  VALUES  FOR  THE 


C I MENSION  FNEMP (3,3) 

OAT A  FNEMP /. *271*00,.* 190 3D0,. 40 816DQ, 

1  .44865C0, .4355300,  .4211600, .486 59D0 , 

2  .*657900, .4448700/ 

CQMMQN/NCRPRN/MEAN ,VARI ,SVARI 
EXTERNAL  NORM , UN  I F 
T  =  1 

100  PRINT  1 

READ  2,MEAN,VARI 
IFIMEAN+VARI.LE.O.CO )STOP 

2  FORMAT ( 2E 5 .0 ) 

SVARI  =  CS CRT ( VAR  I  ) 

WRITE  (4, 14) 

14  FORMAT ( 1H1  ) 

1  FORMAT !  '  +1 NPUT  MEAN  AND  VARIANCE,  FORMAT ( 2F5.0 )• ) 

PRINT  17 

READ  12»WT1,MEAN1 , VAR1, MEAN2 ,  VAR  2 
WT2  =  1  CO  ~  WT1 

17  FORMAT! »+  INPUT  WT1  ,  M5AN1  ♦  VA  R1  ,MFAN2  ,VAR2  ,  POP.MAT  <5  F5  .0  )  •  ) 

12  FORMAT!  5  F5.0) 

PRINT  3 

READ  5, ISEED, NR 
ISEED1  =  ISEED 
DSEED  =  CFLOAT(ISEEC) 

5  FORMAT!  110,15) 

3  FORMAT! '+INPUT  SEEC  AND  NUMBER  OF  REPLICATIONS,  FORMAT! 110,1 5) ' ) 
CO  205  10*1,3 

DO  205  N 1=  1,  3 
NCHI ( N1 , I C )  =  0 
NCHIN!N1,I0)  =  0 
NKS1N1, IQ)  =  0 
NKSN ! N1 , IQ )  =  0 
NFN! N1 , IQ)  =  0 

FNTST ! N1 , I C )  =  ZALC 1 *ZALF ( IQ ) /OSORT! DFLOAT! NSMP! N1 )+ 1 ) )  +  ZALC2 
FNTST !  N1 ,1  C)  =  FNEMP  !M,IQ) 

205  NFNN! N 1 , IQ  )  =  0 
DO  300  1=1, NR 
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CALL  M I XNRM( 50,RN,RD  ,DSSFD,WT1  ,MEAN1 , VAR1 , ME AN2 , VAR2 ) 

CO  250  M  =  l»3 
N2  =  NST(Nl) 

N3  =  NSMP(Nl) 

CALL  VSRTAD(RD(N2),N3) 

CALL  VSRTAD(RN(N2),N3) 

CO  210  NSR=1»N  3 
RNN(NSR) =  RN(NSR  -»  N2-1) 

210  R( NSR l  =  RD( NSR+N2-1 ) 

IDF  =  0 

CALL  GFITCUNIF  ,NCEL  (  N1 )  ,RNN,  N3  ,CELLS  , COMP ,  CHI ,  IDF,Q,  I5R) 

CO  241  10=1,3 

IF  (Q.LT  .CTST  ( IC  )  INCH  I  (Nl,  IQ)  ..=  NCHI(N1»I0)  +  1 

241  CONTINUE 

CALL  NKS1(UNIF,RNN,N3,PDIF,!ER) 

DO  242  IC  *  1.3 

IF(P0IF(1)  .GT.STST(M,IQ)  )NKS(M,IC)  =  NKS(N1,IG)  +  1 

242  CONTINUE 

CALL  FGUTZ(UNIF ,PNN,N3,Q) 

DO  243  10=1,3 

IF(O.GT.FNTST(M,  IC)  >NFN(N1,  IQ)  =  NFN(Nl.IQ) 

243  CONTINUE 
IDF  =  0 

CALL  GF IT ( NORM, NCEL ( Nl ) ,R , N3 , CELLS ,COMP ,CHI , I DF* Q, IER ) 

DO  244  IC  =  1,3 

IF(Q.LT.QTSTUQ)  )NCHIN(N1,IQ)  =  NCHIN(N1,IQ)  +  1 

244  CONTINUE 

CALL  NKS1(N0RM,R  ,N3  .PDIF.ISR) 

CO  245  10=1,3 

IF(PDIF(1).GT.STST(M,IQ))NKSN(N1,IQ)  =  MKSN(N1,I0) 

245  CONTINUE 
CALL  FOUTZ (NCR  M, R , N3  » Q  ) 

DO  246  10=1,3 

IF ( 0 ,GT  *  PNTST( Nl ,  I C  ) ) NFNN( Nl » I  0 )  =  NFNN(N1,IQ»  +  1 

246  CONTINUE 
250  CONTINUE 
300  CONTINUE 

ISEED  =  CSEED 

WRITE  (  4, 7)  ME  AN,  VAR  I  ,  NR,  I  SEED  1  ,ISEEC 
WRITE (4, 11  )WT1,MEAN1,VAR1,WT2,MEAN2,VAR2 
11  FORMAT ( //'  MIXED  NORMAL  (MIXNCPM)  =  • , F5 .2 , • *N ( • , F3 . 1 , • , • , F 3 . 1 , 

1  ')  +' *  F5. 2,  **N(*,F3.1,*,*,F2.1,*  )•//) 
WPITE(4,8)CHISCU,((NCHI(N1,!Q). 

1  Nl=l, 31  ,10  =  1,31  ,((NCHIN(N1,IC), 

2  Nl=l,  3  ) ,  10=1,  3  ) 

WRITE(4,8)K0LSMR,( (NKS(N1,IQ) ,N1=1,3), 10=1,3 », 

1  (  (NKSN(M.IQ)  ,N1  =  1  ,3)  ,10=1,3  ) 

WRITE(4,8IF0UTZL, ( (NFN(N1, IQ) ,N1=1,3) ,10=1,3) , ( ( NFMN ( N1 ,1 0 ) ,M1=1 
1  )  ,  10  =  1  ,3) 

7  FORMAT ( // '  THE  VALUES  OF  MEAN  (U)  AND  VARIANCE  (V)  ARE* ,2F12.5// 
9  •  THE  RESULTS  OF  *,I5,'  REPL !CAT ICN S,  STARTING  WITH  SEED', 

1  I12//33X, 'ENDING  WITH  NEXT  SEEC',112//) 

8  FORMAT  (///20X,  2A6//2CX,  '20  PTS*, 

2  6X , ' 30  PTS'  ,6  X , ' 50  PTS • //10X ,  •  10  %  • , 31 12/ •  CONTROL  % 2X, •  5T*, 

3  31 12/10X, '  1?* ,2112///'  MIXNCPM* ,2X, '10** ,3112/'  VS',5X,»  5*' 

4  3112/'  N  (U  ,V  ) '  ,  2X ,  '  1?  *,  31 12  ) 

WRITE (4, 9) (( FNTST(N1 ,10) ,N1 =1 , 3 ) , I C=1 ,3 ) 

9  FORMAT ( / *0  THIS  WAS  BASEC  ON  THE  FOLLLOWING  TEST  VALUES*, 

1  / (13X , 3F12.5 I  ) 

GO  TO  ICO 
END 


SUBROUTINE  M IXNRM (NO , RAN,R ,DSEED ,WT1 ,MEAM1 , VAR  1 ,MEAN2 , VAR 2 ) 
IMPLICIT  REAL*8  (A-M,0-Z> 

DIMENSION  R ( 1 )  ,RAN(1) 

REAL*4  RN(IOO) ,  RNN  » RNNN 
SVAR1  =  DSQRT(VAFl) 

SVAR2  =  DSQRT ( VAR  2 ) 

WT2  =  l.CO  -  WT1 
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200 


XX  =  X*  X 

TAN2  =  (P0  +  XX*(P1  ♦  XX*P  2) ) /(PIBY4*(QQ  +  XXMQ1  +  XX*02))» 
RETURN 

TAN2  =  TAN(XARG) / XA  PG 

RETURN 

END 


200 

300 

400 


SUBROUTINE  FCUTZ  (FCCF,XT»NXT»FN) 

THIS  SUBROUTINE  GENERATES  THE  STATISTIC  FOR  THE  FCUTZ 
INPUT  VARIABLES  ARE 


FN  TEST, 


PCDF  -  THE  CUMULATIVE  DISTRIBUTION  FUNCTION  AGAINST  WHICH  THE 

DEVIATES  ARE  BEING  TESTED.  CALLING  SEQUENCE  MUST  BE  OF 
THE  FCRM  •  CALL  PCDF(X,P)  ',  WHERE  X  IS  AN  INPUT  VALUE, 
AND  THE  VALUE  OF  THE  CUMULATIVE  DISTRIBUTION  FUNCTION 
IS  RETURNED  IN  P. 

P  MUST  EE  BETWEEN  0  AND  1. 

XT  -  THE  ARRAY  OF  DEVIATES,  IN  INCREASING  ORDER. 

NXT  -  THE  NUMBER  OF  DEVIATES  IN  TH5  ARRAY  XT  <  *  N  -  1 > 

THE  RETURNED  VALLE  IS  FN,  THE  VALUE  OF  T HE  STATISTIC. 


NXT  I 


„  ,r  =NTLY  LIMITED 

THE  VARIABLE  XTD. 


EB£S_  x 

CIMENS  ION  XT ( 1 ) 

REAL*8  XTD ( 51 ) ,P N , FN D 
N  =  NXT  ♦  1 
DO  200  1  =  1,  NXT 
K  =  N  -  I 

CALL  PCDF(XT(K)»P) 

XTD ( K+l )  =  P 

RN  =  l.DC/N 

XTD(l)  =  RN  -  XTD ( 2 ) 

DO  300  I =2, NXT 

XTD(I)  =  RN  -  XTC ( I +1 )  +  XTD ( I ) 
XTD(N)  =  RN  -  l.CC  +  XTD( N ) 

FND  =  0. 

CO  400  I  =  1  ,N 

FNO  =  FNC  ♦  CMAX1(XTC(I),0.00) 

FN  =  FND 

RETURN 

END 


TO  A  MAXIMUM  OF  50  BY  THE  DIMENSION  OF 
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150 

200 


CALL  GGUBS(DSEED,2*ND,RN) 

CO  200  NN=1 , NO 
RAN(NN)  =  PN(NN) 

CALL  MDNRI$(RN(NN),RNN,IER) 
IF(RN(NN4ND).GT.WT2 )G0  t0  150 
R(NN)  =  PNN*SVAR2  +  MEAN2 
GO  TO  200 

R ( NN )  =  RNN+SVAR1  +  MEAN1 

CONTINUE 

RETURN 

END 


SUBROUTINE  NCRR(XtP) 

PEAL*8  MSAN.VARI ,SVARI 
COMMON/ NOR PRM/ME AN  »VARI,SVAPT 
T  =  (X  -  MEANI/SVAPI 
P  =  .5*£RFC<-T*. 1071068) 
RETURN 
END 


SUBROUTINE  UNI  F  <  X  ,  P I 

P  *  X 

RETURN 

END 


200 

300 

400 


TS 


EST. 


SUBROUTINE  FOUTZ  ( PCOF,XT,NXT , FN ) 

THIS  SUBROUTINE  GENERATES  THE  STATISTIC  FOR  THE  FCUTZ  FN 
INPUT  VAPIABLES  APE: 

PCD F  -  THE  CUMULATIVE  DISTRIBUTION  FUNCTION  AGAINST  WHICH  THE 
DEVIATES  ARE  BEING  TESTED.  CALLING  SEQUENCE  MUST  RE 
THE  FORM  •  CALL  PCCF(X,P|  »,  WHERE  X  IS  AM  INPUT  VALU 
AND  ThE  VALUE  OF  THE  CUMULATIVE  DISTRIBUTION  FUNCTION 
IS  RETURNED  IN  P. 

P  MUST  eE  EETWECN  0  AND  1 . 

XT  -  THE  ARRAY  OF  DEVIATES,  IN  INCREASING  ORDER. 

NXT  -  THE  NUMBER  OF  DEVIATES  IN  THE  ARRAY  XT  ( *  N  -  1) 


THE  RETURNED  VALLE  IS  FN,  TH£  VALUE  OF  THE  STATISTIC. 


NXT  IS  PRESENTLY  LIMITED  TO  A  MAXIMUM  OF  50  BY  THE  DIMENSION  OF 
THE  VARIABLE  XTD . 

DIMENSION  XT  (1  ) 

REAL*8  XTD ( 51 ) ,RN ,FNO 
N  =  NXT  4  1 
DO  200  1=1, NXT 
K  =  N  -  I 

CALL  PCCF(XT(K),P) 

XTD(K+1>  =  P 

RN  =  l.DO/N 

XT D ( 1 )  =  RN  -  XTC ( 2 ) 

DO  300  I =2, NXT 

XTD(I)  =  RN  -  XTO ( I +1 )  4  XTD ( I ) 

XTD(N)  =  RN  -  l.CO  4  XTC(N) 

FND  =  0. 

CO  400  1=1, N 

FND  =  FNC  4  CMAX1(XTD(I),0.D0) 

FN  =  FND 

RETURN 

END 
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THIS  PROGRAM  GENERATE S  P? ARSON  TYPE  I  OR  II  RANDOM  DEVIATES. 

THE  PARAMETERS  ARE  CALCULATED  IN  TERMS  OF  81  AND  B2  IN 
SUBROUTINE  PRM.  SUBROUTINE  ACINT1  IS  USKD  tq  CALCULATE  THE 
CDF  TABLE.  GGUBS  IS  USED  TO  GENERATE  RANDOM  (O.ll  DEVIATES 
(UNIFORMLY  DISTRIBUTED,  AMD  THEN  SUBROUTINE  RANDP1  DOES 
AND  INVERSE  CALCLLATICN  TO  OBTAIN  THE  RANDOM 
PEARSON  CEVIATE.) 

A  SET  OF  50  RANDOM  CEVIATES  ARE  GENERATED,  THEN  THE  CHI 
SQUARED  TEST,  THE  KCLMOLGOROV-SMI RNOV  TEST,  AMD  THE  FOUTZ  <=N 
TEST  ARE  APPLIED.  THE  UNIFORM  DEVIATES  ARE  TESTED  AGAINST 
THE  HYPOTHESIS  THEY  CAME  FROM  A  UNIFORM  DISTRIBUTION  AS  A  CONTROL 
THEN  THE  PEARSON  OEVIATES  ARE  TESTED  AGAINST  THE  HYPOTHESIS 
THAT  THFY  CAME  FPCM  A  NORMAL  DISTRIBUTION. 

THIS  TEST  IS  REFLICATEO  A  NUMBER  OF  TIMES. 

THE  TESTS  ARE  APPLIED  TO  THE  FIRST  20,  tHc  NEXT  30,  AND  ALL  50. 


THE  SET  CF  DEVIATES  IS  ALSO  TESTED  AGAINST  THE 
THIS  IS  REPEATED  FOR  THE  SAME  SETS. 


NORMAL  DISTRIBUTION 


THE  INPUT  VALUES  ARE  B1.B2, I $EC0,  AND  THE  NUMBER  OF  REPLICATIONS 
IMPLICIT  REAL*8  (A-Z) 

INTEGERS  T,IER,N,  I ,  NP1,N  1 ,  N  2  ,N3 , 1  SEED ,  NR  ,NSTOP ,  I  DF  ,  I  Q  ,  I  SEED! 

1  ♦ NCEL ( 3 ) , NST ( 3 ) , N$MP (3 ) , NCH I ( 3 ,3 ) , NCHIN ( 3, 3 > , NKS ( 3, 3 > , 

2  NKSN( 3 » 3 1 »NFNN( 3 ,3)  ,NFN(3,3) , NAFNN ( 3 ,3 ) , NA'N ( 3 , 3 ) , NS R 

DI MENS  ION  X( 2001  ) ,COF( 2001 ) , A ( 2001 ) , R ( 2001 » , RD ( 55) ,CHISQU(2) , 

1  KOLSMR  ( 2 )  »P0UTZL(2)  ,ASFL(2)  ,  f=NTST  (3 ,3  ) ,  ZALF  (  3  ) ,  RN(  5  5  ) 

COMMON  /PPARM/C0,C1,C2, Al, A2,BCC,BC1,M1 ,M2,K!NV,XL,XR,MEAN,T 
COMMON  /CDF/X,CCF,A,B,NP1 

REAL*4  CELLS (50)  ,CCM P ( 50) , Q , CHI , P ( 55 ) , PD  I F (6 ) , QTST ( 3 ) , RNN ( 55 ) , 

1  STST (3,2) 

DATA  STST/. 26473,. 21 756,. 16959,. 29408, .24170, . 18841, . 35241, 

1  .28987, .22604/ 

CATA  NCEL/ 8, 10 , 10/, NST/1, 21 , 1 / ,NSMP/ 20, 30 , 50/ , QTST/. 1 ,. 05 , . 01 / 

1  ,CHISCU/'CHI  SQUA' , ‘RED  TEST*/. KCLSMR/ »KOLMOG-S • , 'MIR  TEST'/ 

2  , FOUTZ L / ' EMP  FOLT'.'Z  TEST* /, ZA LC 1 , ZALC2/ .243 06  9C0 , .367879  00/ 

3  , ZALF/1.28155C0, 1.6  448  5D0,  2. 3 26 3 SCO/ 

DIMENSION  FNEMP (3 , 3  ) 

DATA  FNEMP/. 42714D0, .4190300, .4081600  , 

2  .4436500, .4355300, .4211600  , 

3  .4865900, .4657900, .4448700/ 

4  , AS FL/ ' AS YM  FOU'.'TZ  TEST’/ 

EXTERNAL  NORM , UN  I F 

T  =  1 

100  PRINT  1 

READ  2,B1,B2 
IF(B1+B2.LE.0.D0) STCP 

2  FORMAT (2E5.0) 

WRITE (4,14) 

14  FORMAT ( 1 H 1 1 

1  F0RMAT(»4  B1  AND  B2,  FORMAT ( 2= 5.0 )• ) 

CALL  PRM(B1,82,1.C0,IER) 

IF ( IER . NE . 0 ) GO  TG  SCO 
A(l)  =  0.00 
CDF(l)  =  0. 

CALL  INTSZ(N.X) 

NP1  =  N  ♦  1 
PRINT  3 

READ  5,  ISSED.NR 
ISEED1  =  I  SEED 
DSEED  =  CFLOATd  SEED) 

5  FORMAT (110, 15) 

3  FORMAT ( * +1 NPUT  SEED  AND  NUMBER  OF  REPLICATIONS,  FORMAT( 110, 15 ) • ) 
CO  205  I  Q=  1, 3 

DO  205  Nl=l , 3 
NCHI (Ml ,IQ)  =  0 
NCHIN(NI.IQ)  =  0 
NKS (Ml, I C)  =  0 
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205 


210 


241 


242 


243 


244 


245 


246 

250 

300 


+  ZALC2 


NKSN(NltlO)  =  0 
NFN  ( N1 1  10=0 

FNTST< N1 *10)  =  ZALC1*ZALF(10) /DS0RT(DFL0AT(NSMP(N1)+1) ) 

NAFNN<  N1,IQ)  =  0 
NAFN(N1,IQ)  =  0 
NFNN<N1,IQ)  =  0 
DO  300  I =1 »NR 

CALL  RANPD1(50»RN,RC,0SSFD) 

DO  250  Nl=l»  3 
N2  =  NST(Nl) 

N3  =  NSMP(Nl) 

CALL  VSPTAD<RD(N2) ,N2) 

CALL  VSRTAD(RN(N2 ),N2) 

DO  210  NSR=1 »N3 
RNN(NSR)=  RN ( NSR  ♦  N2-1) 

R( NSR I  =  RD(NSR+N2-1> 

IDF  =  0 

CALL  GFIT(UNIF,NCEL(N1),RNN,N3,CELLS,C0MP,CHI ,  IDF  tQt I ER) 

CC  241  10=1,3 

IF(0,L7.0TST<IQ) INCHI(N1,IQ)  =  NCHI(N1,IQ)  +  1 
CONTINUE 

CALL  NKSKUNIF  ,  RAN,  N3,  PDIF,  IER) 

DO  242  10  =  1,2 

IF(P0IF(1I  .GT.STST(Nl,IO)NKS(Nl,  IQ)  =  NKS(N1,IQ)  +  1 
CONTINUE 

CALL  F0UTZ(UNIF»RNN»N3»Q) 

DO  243  I C=1 r 3 

IF<Q.GT.FNEMP<N1,IC)  )NFN(N1,IQ>  =  NFN(N1,I0>  +  1 
I F (0, GT • FNTST ( N1 ,  10  )  ) NApN ( N1 , IQ )  =  NAFN<  N1 » I Q)  +  1 
CONTINUE 
IDF  =  0 

CALL  GFIT(N0RM,NCEL(N1),R,N3, CELLS, COMP, CHI , I  DP, Q, IER ) 

DO  244  IQ  =  1,3 

IFIO.LT  .OTST(IQ) )NCHIN(N1,!Q)  =  NCHIN(  N1 ,10)  +  1 
CONTINUE 

CALL  NKS1(N0RM,R,N3, PDIF, IER) 

DO  245  10=1,3 

IF(PDIF(1).GT.STST(M,IQ))NKSN(N1, IQ)  =  NKSN(N1,IQI  ♦  1 
CONTINUE 

CALL  F0UTZ(N0RM,R,N2,Q) 

DO  246  10=1,3 

IF(0.GT.FNEMP(N1,IC))NFNN(N1,!C)  =  NFNN(N1,I0)  +  1 
IF(0.GT*FNTST(N1,IC) )NAPNN ( Ml , 10 )  =  NAFNN( N1 , IQ)  +  1 
CONTINUE 
CONTINUE 
CONTINUE 
ISEED  =  DSEED 
WRITE (4 ♦ 7  )  B1 , B2, A  1, A  2, Ml, M2, NR, I  SEED  1,1  SEED 
WRITE(4,8ICHIS0U,<  (NCHI  (N1 ,10)  , 

1  N 1=1, 3  )  ,10=1, 3 ) , ( (NCHIN(NI.IQ), 

2  N 1 = 1 , 3 ) ,10=1,3) 

WRITE(4,£)KCLSf'R,(  (NKS(N1,I0)  , Nl=l  ,3  ) ,  10=1 ,3  )  , 

1  ( <NKSN(N1,IQ) ,N1=1,3),IQ=1,3) 

WRITE (4,81 FOUTZL, ( (NFN(N1, IQ) ,N1=1,3), 10=1,3), ( ( NFNN ( N1,IQ ) ,N1=] 
1  ),  10=1,3) 

WR ITE (4 , 8 ) AS  FL , ( (  NAFN <N1, 1 0 ) ,  Nl=  1,  2) ,  IQ  =  1,  3) , 

1  ( (NAFNN(N1,IC) , N1 =1 ,3 ) , 1 0=1 ,3 ) 

7  FORMAT ( / / •  THE  VALUES  OF  B1  AND  B2  ARE • ,2F12 ,5 /// 

A  •  VALUES  0=  Al,  A2,  Ml,  AND  M2  ARE* //4F12. 6/// 

9  •  THE  RESULTS  CF  ',15,*  REPLICATIONS,  STARTING  WITH  SEED*, 

1  I12//33X, ‘ENDING  WITH  N^XT  SEED*, 112//) 

8  FORMAT (///20X,2A8//20X,‘20  PTS‘- 


,6X,‘5C  PTS'//10X,‘10%‘ ,3112/' 

1** ,2112///'  PEARSON  ' , 2X , • IOX',3112/* 
15?'  ,3112) 


2  6X,'30  PTS 

3  3I12/10X, • 

4  3112/' 

GO  TO  100 

900  PRINT  4,N,IFR,B1,E2,ERR0R 

4  FORMAT ( ' N  »I5R, 61 »B2 '  ,2 110 , 1 P3 C15 .6 ) 
STOP 
END 


5?' 


VS' ,5X, '  5* ' 
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SUBROUTINE  RAN  PCI  (NIC  ♦  RAN.R  t  DSEFD  ) 
IMPLICIT  REAL* 6  (A-M.C-Z) 

CIMENS ION  R  ( 1 ) »R  AN  < 1 ) 

REAL*4  RMlOOItRNN 
CALL  GGUBS(DSEED ,NC,RN) 

DO  200  NN=1»ND 
RAN<  NN I  =  RN(NN) 

CALL  CDF  INV(  RAN(  NN  )  t  RNN) 

R(NN)  =  PNN 
CONTINUE 
RETURN 
END 


SUBROUTINE  NORMX.F) 

P  =  .5*ERPC(-X*.7C71C68) 

RETURN 

END 


SUBROUTINE  UNI F(X,P) 

P  =  X 

RETURN 

END 


SUBROUTINE  PRM ( B1 t E2 tU2 » IER ) 

IMPLICIT  REAL*  8  (A-Z) 

COMMON  /PPARM/C0.C1.C2»A1, A 2 » BIGCO »BI GC1 , Ml ,  M2 t K IN V, XL  ♦  XR  »  MEAN  » T 
INT£GSR*4  T»IER 
IER  =  0 

I  F (T.EO • 2  I  82  =  (6.  +  3.*Bl)/2. 

IF  (  T.E  C  .  5  )  31  =  4.* (2. *82  -  6 . ) * (4 . *B2-3 . I / ( ( 82+3 . )**2  + 

1  12.*<4.*B2  -2.)) 

DEN  =  10.*B2  -  12.*ei  -  18. 

CO  =  ( 4 • *B2  -  3.*ei)*U2/DEN 
Cl  =  DS0PT(U2*ei )*<B2  +  3.)/DEN 
C2  =  (2.*B2  -  3.*E1  -  6.J/DEN 
A1  =  0. 

A2  =  0. 

BIGCO  =  0. 

BIGC1  =  C. 

Ml  =  0. 

M2  =  0. 

GO  TO  (100,200,30C,400,500,600,7CO)tT 
100  CALL  RCCTS (C0»CltC2.Al, 42,  IER) 

IF(IER.NE.O)RETUPN 

Ml  =  (Cl  +  A1)/C2/(A2  -  Al) 

M2  =  —(Cl  +  A2)/C2/(A2  -  AM 

KINV  =  <A2  -  A  1 1 **( M1+  M2  +  1 . )*DGAMMA (Ml+1 . ) *OGANMA ( M2+1 . ) 

1  /DGAMMA (Ml  +  M2  *  2.) 

XL  =  Al 
XP  =  A2 
GO  TO  750 
200  GO  TO  100 
300  Ml  =  (C0/C1  -  Cll/Cl 
C2  =  0. 

KINV  =  DGA  NMA( Nl  +  1 • ) ADEXP (C0/C1**2 )*DABS ( Cl ) **( 2 . *M1  +  1 . ) 

XL  =  -C0/C1 
XR  =  1.C50 
GO  TO  75C 

400  BIGCO  =  CO  -  Cl**2/C2/4. 

BIGC1  =  C1/C2/2. 

Ml  =  (Cl  -  BIGC1  )/0$QRT(2. DO) /BIGCO 
M2  =  DSQRT (BIGC0/C2 ) 

KINV  =  1.00 
XL  =  -1.D50 
XR  =  1.050 
GO  TO  75C 

500  BIGC1  =  C1/2./C2 
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600 


Ml  =  (Cl  -  BIGCD/C2 
IF(C2.GT.1.0C)G0  70  800 

KINV  =  CABStMl )**(1./C2  -  1 .  > /DGAMMA  ( 1.  /C  2-1.  ) 

XL  =  -BTC-C1 
XR  =  1 .050 
GO  TO  750 

CALL  R0CTS(C0*Cl»C2tAltA2»IER) 

Ml  =  -(Cl  ♦  A 1 ) /C 2 /( A2  -  Al) 

M2  *  (Cl  ♦  A2)/C2/(A2  -  Al I 

IF(M2.GE.-1.D0.0R.M1  +  M2 . G5 .0  .00 » GO  TO  800 

KINV  =  ( A2  -  A1)**(M1  +  M2  +  1. ) *D GAMMA ( M 2+1. ) *DGAMMA ( -Ml -M2-1 • ) 

1  /DGAMMA (-Ml) 

XL  =  A  2 
XR.  =  1.050 
GO  TO  750 

700  KINV=C0**(-.5D0/C2)*DSQRT(C0/C2)*DGAMMA(.5D0)*DGAMMA<  .5D0/C2-.5C 
1  /D GAMMA (.5D0/C2  ) 

Cl  =  0. 

XL  =  -1.050 
XP  =  l  C50 

MEAN  = * ( 4. D0*B 2  -  3 .D0*B1 ) /3. DO/ ( B2  -  B1  -  1 . DO ) *U2*PCF(0 .00 ) 
RETURN 
IER  =  1 
RETURN 
END 


750 

800 


500 


600 


SUBROUTINE  ROOTS ( CO, C1,C2, A 1 , A 2, IER) 
IMPLICIT  REAL* 8  (A-Z) 

INTEGER*4  IER 
IER  *  0 

IF(C2.EC.0.D0)GC  TO  500 

DIS  =  C 1**2  -  4. CC*C0*C2 

IF(DIS.LT.0.D0»GC  TC  600 

SDIS  =  0  SORT ( D I S I 

DNUM  =  -Cl  -  SDIS 

IF(C1.LT.0.D0)CNUM  =  -Cl  +  SCIS 

X2  =  DNUM/2.DO/C2 

XI  =  C0/C2/X2 

Al  =  DMIM(X1,X2) 

A2  =  DMAX1 ( Xl»  X2 ) 

RETURN 

A2  =  1.D75 

Al  =  -C0/C1 

RETURN 

IER  =  -1 

RETURN 

END 


FUNCTION  PDF (X ) 


THIS  FUNCTION  EVALUATESTHE  PEARSON  DISTRIBUTION 
FOR  A  GIVEN  X,  THE  PARAMETERS  HAVING  PREVIOUSLY  8EEN 
CALCULATED  IN  SLBRCLTINE  PRM. 


100 


IMPLICIT  REAL*8  (A-2) 

INTEGER*^  ITY 

COMMON/ PPARM/  CO t Cl t C2» A 1 , A 2 » B IGCO t BI GC 1 , M 1 »M2 »K I NVt XI , X2 t MEA M « I 
GO  TO  (lQOt200»300»400t500»600t70Q)t!TY 
IF( X.LE.A1.0R.X.GE.A2)GO  TO  140 
PDF  =  (X  -  A1I**M1*(A2  -  X )**M2/K INV 
RETURN 

140  IF(X.LE.A1.AND.M1.LT.0.D0»G0  TO  150 
IF(X.GE.A2.AND.M2.LT.0.C0)G0  TO  15C 
PDF  =  1.0-25 
RETURN 

150  PDF  =  1.C25 
RETURN 

200  C-0  TO  IOC 
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300  PDF  =  (CO  +  C1*X  )**I*1*DEXP(-X/C1  )/KlNV 
RETURN 

400  PDF  =  (BIGCO  +  C 2* ( X+BIGC 1 ) **2 ) ** ( -1. /C 2) * 

1  DEXP  ( -M*DATA  N  (  (X  +  P IGC1 ) / M2 ) ) / K  INV 
RFTURN 

500  PDF  =  (X  +  BIGC1  )**(-l./C2  )*DSXP(Ml/(  X  ♦  BIGCDI/KINV 
RETURN 

600  PDF  =  (X  -  A 1 ) 1* ( X  -  A2)**M2/KINV 
RETURN 

700  PDF  =  (CC  ♦  C2*X**2)**(-.5D0/C2)/KINV 
RETURN 
END 


150 

160 


170 


300 


301 


302 


304 


308 


309 

310 


400 


SUBROUTINE  INTSZ(N.X) 

IMPLICIT  REAl*8  (A-Z) 

COMMON  /PPARM/C0»C1,C2, Al?A2t 8C0tBCllMl,M2tKINVt XL,XR,M=ANtT 
COMMON /CCF/XDUMC 2001) , CDF ( 2001 ), A ( 2001 ), B ( 2001  ), NP1 
Cl  MENS ICN  X(l) 

INTSGER*4  T.N, NN , I  NT , IEND » NP1 ,  IER 
X(l)  =  A 1 
SCI)  =  PCF(Al) 

CDF ( 1 )  =  O.DO 
INT  =  1 
NN  »  1 

DXM  =  (A2  -  A1J/1.C7 
CX  =  DXM*1  *03 
GO  TO  160 

DX  =  D  X* ( 1 .DO  +DSCRT( .475D-4/-REST ) ) /2. 

INT  =  1 

XCNN+l )  =  X(NN)  4  CX 
IEND  =  1 

IF(X(NN  +  1) .LT.A2  )G0  TO  170 

X( NN+l )  =  A2 

DX  =  X ( NN+l)  -  X ( N N ) 

IEND  =  2 

B(  NN4-1 )  =  PDF(  X(  NN4-1 )  ) 

ESTINT  =  ( B ( NN  4 1 )  4  B (NN ) ) *D X/ 2. DO 

IF(ESTINT.LT.l .C-8 )GQ  TO  301 

DC DF  =  ADINTK  X(NN)  ,X(NN4-1)  tFRRCP.»ISR) 

IF ( IER. GT. 100) GO  TO  900 
IF(DCDF.GT.,025D0)GC  TO  309 

EREST  =  (l.D0/B(NN4l)  -  l.DO/8 (NN) )*DCDF/8 .DO 
A( NN)  =  EREST 
EREST  =  CABS(EREST) 

GO  TO  (  2C2.308),  INT 

DCDF  =  ESTINT 

EREST  =  5. 27777777 77D-6 

A ( NN )  =  EREST 

GO  TO  308 

IF(ERE$T.LE..5C-4)G0  TO  308 
CX  =  DX*(DSQRT(.25C-4/EREST) 

IF(DX.GT.DXM)GC  TC  304 
DX  =  DXM 

X(NN4-1)  =  X  ( NN )  4  CX 

INT  =  2 
GO  TO  17C 
X(  NN4-1  )  = 

IEND  =  1 
GO  TO  170 
CDF(NN4-1) 

NN  =  NN  4 


4  l.DO*.  5D0 


X ( NN )  4CX 


4  DCDF 


910 


=  CDF(NN) 

1 

IF (NN. GT .2000 ) GO  TO 
GO  TO  ( 150.310) .IEND 
DX  =  DX*.0125DG/CC0F 
GO  TO  160 
N  =  NN  -  1 
NP1  =  NN 
DO  400  NN=1 » N 

CDF(  NN4-1 )  =  CDF(NN4-1)/C0F(N4-1) 
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2  FORMAT ( '  CDF(NN)  =  *,F20.15) 

RETURN 

900  PRINT  1,NN»X(NN) ,X(NN+1) »EFRCR 
STOP 

1  FORMAT ( *  I »XI»XI+1,ERPGR,»I5,1P3C15.6) 
910  PRINT  4 , X( NN  ) 

4  FORMAT ( '  RUN  CUT  OF  SPACE  AT  X=',F12.5) 
STOP 
ENO 


3  00 


350 

360 

400 

1 

410 


SUBROUTINE  CDFINV(C,P) 

REAL*8  X,CDF,A, B , CC , DCC,A0INT1 ,  XC, E 
COMMON/CDF /X( 2001)  ,CCF (2001 )  , A (2001 ), B (2001 ), NP1 
CALL  TABL0C(CDF»CELE(C),NP1,I  ) 

XC  =  X(I)  +  (C-CCFU  )  1/ (CDF(  1  +  1)  -  CDF(I))*(X(  1  +  1)  - 
IF(DABS(A(I)).LE..5D-4)G0  TO  410 
CO  400  J=l,10 

OC  =  AOINTKX(I)  ,XC,E,IEP) 

OOC  =  ( CDF ( I )  ♦  DC  -  C)/PDF(XC) 

XC  =  XC  -  DDC 

IF (DAB  S(CDC) .LT..5C-4)G0  TO  410 
IF  (  XC  .  LE  .X  ( I  )  )  GO  TO  350 
IF  (XC.GE.X(I+1  )  )GC  TO  360 
GO  TO  4CC 

XC  =  X(I1  +  (C  -  COF(I) )/DC*< XC  -  XC I » > 


X(  !  >  ) 


SP  J°xWi 


X(  1+1) 

CONTINUE 
URI TE  (  6  » 1 )  C  ,  XC 
FORMAT ( •  10  ITERATIONS 
p  -  xf 
RETURN 
END 


+  (C  -  CDF(I+1) ) /(COF(I)  +  DC  -  CDF( 1+1 ) ) * (X ( 1+1 )-X 


IN  CDF INV ,  C  AND  XC  =«,2E12.4> 


SUBROUTINE  TABLOC( XT, X, M,NL ) 
REAL*8  XT ( 1 ) ,  X 

NT  =  ALGG(FLOAT(M  ) /.301  +  1. 
NU  =  M 
NL  =  1 

DO  200  1=1, NT 
NG  =  (NU  +  NL)/2 
IF ( X.GE • XT (NG)  I G C  TC  100 
NU  =  NG 
GO  TO  200 
100  NL  =  NG 
200  CONTINUE 
RETURN 
END 


FUNCTION  ADINTKXl  ,X2  ,ERROR,  IER) 

THIS  SUBROUTINE  USES  DCADR5  TO  OBTAIN  THE  INTEGRAL  OF  A 
PEARSON  TYPE  I  OF  II  CISTR I5UT ION  FUNCTION.  SUBTRACTING 
OUT  THE  SINGULARITY  IS  USED. 

IMPLICIT  REAL*8  (A-Z) 

INTEGER*4  T,  IER 
EXTERNAL  F1,F2,PCF 

COMMON /PPARM/CO ,C1»C2»A1»A2,BC0,8C1»M1,M2»KINV»XL ,XR , MEAN ,T 
IF ( X 1 .GT  .(XL+XR)/2.CQ)G0  TO  2C0 
IF ( Ml . GT .0 .DO ) GC  TC  300 

PR  =  DCADRE(F1, XI, X2,l.C-6, 0.00, ERROR, IER) 

ADINT1  =  PR  +  <A2-A1)**M2*(  (  X2-A1 )  **  <  M1+ 1.  DO)  -  (  Xl-Al )  **  ( MH-1.  DO 
1  /(Ml  +  l.CO)/KINV 
RETURN 

200  I F ( M2 . GT  .0 .DO ) GO  TC  300 

PR  =  DCACRE(F2,X1,X2,1. D-6 ,O.CO,ERFOR, IER) 

ADINT1  =  PR- (A 2  -  A1)**M1*((A2-X2)**(M2+1.00)-(A2-X1)**(M2+1.DO 
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1  / (M2+1.C0 )/KINV 
RETURN 

300  ADINT1  =  0CADRE( P0F,X1,X2, 1 . D-6, 0. CO ♦ ERROR , IEP ) 
RETURN 
END 


FUNCTION  FKX) 

IMPLICIT  RSAL*fi  (A-Z) 

COMMON  /FPARM/C0,C1,C2,A1,A2, B CO  ,  BC1 ,  Ml , M2, K  IN  V ♦  XL , XR,MEAN,T 
IF ( X.LE  «A1 )G0  TO  2  0  C 

FI  =  T  X-Al  )**M1*  (  (A2  -  X l **M2  -  (A2  -  Al> **M2) /KINV 
RETURN 

200  FI  =  O.DC 
RETURN 
END 


200 


FUNCTION  F2(  X) 

IMPLICIT  REAL*8( A-Z ) 

CCMMON/P FARM/CO t Cl  , C2 , A1 , A2 , 8C0, 8C 1, Ml , M2, K INV , XL , XR  ,  MEAN  ,  T 
IF ( X.GE  «A2  IGO  TO  2C0 

F2  =  (A2-  X )**M2*( (X-A1)**M1  -  < A 2-A 1 ) **M 1 ) /K I  NV 

RETURN 

F2  =  O.DC 

RETURN 

END 


200 

300 

400 


SUBROUTINE  FCUTZ  (FCOF,XT,NX?,FN) 

THIS  SUBROUTINE  GENERATES  THE  STATISTIC  for  THS  FCUTZ  FN  TEST. 
INPUT  VARIABLES  ARE 


PCDF  - 


XT 

NXT 


THE  CUMULATIVE  DISTRIBUTION  FUNCTION  AGAINST  WHICH  TH? 
DEVIATES  ARE  BEING  TESTED.  CALLING  SEQUENCE  MUST  Bc  OF 
THE  FCRM  •  CALL  PCCF(X,P>  ',  WHERE  X  IS  AN  INPUT  VALUE, 
AND  THE  VALUE  OF  THE  CUMULATIVE  DISTRIBUTION  FUNCTION 
IS  RETURNED  IN  P. 

P  MUST  EE  BETWEEN  0  AND  1. 

THE  ARRAV  OF  DCVIATE$,  IN  INCREASING  ORDER. 

THE  NUMBER  OF  DEVIATES  IN  THE  ARRAY  XT  C*  N  -  1) 


THE  RETURNED  VALUE  IS  FN,  THE 

A 


NXT  IS  PRESENTLY  LIMITED  TO 
THE  VARIABLE  XTO. 

DIMENSION  XT  (1  ) 

REAL*8  XTD< 51) ,RN ,FND 
N  =  NXT  +  I 
00  200  I =1  ,NXT 
K  *  N  -  I 

CALL  PCDF(XT(K),P) 

XTD( K+l )  =  P 

RN  =  l.DC/N 

XTD(I)  *  RN  -  XTCC2) 

DO  300  I *2, NXT 

XTD(I)  *  RN  -  XTD ( 1+1)  +  XTD ( I ) 
XT  D  <  N )  *  RN  -  l.CO  4  XTO(N) 

FND  =  0. 

CO  400  1*1, N 

FND  =  FND  +  CMAX1(XTC(I),0.D3I 

FN  =  FND 

RETURN 

END 


VALUE  OF  THE  STATISTIC. 

MAXIMUM  OF  50  BY  THE  DIMENSION  OF 
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APPENDIX  2 


CHI  SQUAREC  TEST 

20  PTS  30  DTS  50  PTS 


1C3 

410 

4  29 

461 

CONTROL 

5 3 

2  40 

240 

229 

1? 

48 

55 

56 

PS  SF 

10? 

1784 

2908 

4125 

VS 

53 

1362 

2384 

3645 

N  (0,1) 

1? 

532 

1410 

2673 

KOLMOG-SMIR 

TEST 

2  C  PTS 

30  PTS 

5  C  PT 

103 

491 

495 

4  63 

CCNTFOL 

53 

253 

253 

224 

13 

49 

58 

46 

PSS  f 

1C3 

1636 

220  2 

3316 

V  S 

53 

1008 

1401 

2410 

N  (0,1) 

13 

270 

466 

9C5 

FQUTZ  F N  TEST 

2 C  PTS  30  PTS  50  PTS 


103 

525 

481 

488 

CCNTRCL 

5? 

225 

250 

253 

1? 

4C 

35 

62 

RSSF 

10? 

36  03 

4010 

45 9  T 

VS 

53 

2999 

3515 

4279 

N  (0,1) 

1? 

1894 

2463 

346  6 

’CCNTPCL'  IS  TFE  TEST  CF  INI  FORM  (0*1)  VS  UNIFCRM(C,1) 

•  R  S  SF  '  IS  THE  FANDOM  STAEILIZED  STANDARD  FORM 
ALPHA  =  1.00,  6ETA  =  C.O 


NUM8ER  OF  H  VP  C  THE  SI  S  REJECTIONS  IN  5000  REPLICATIONS 
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CCMFCl 


RSS  F 
VS 

N  (  C  ,1  ) 


CONTROL 


RSSF 

VS 

N  (0,1) 


CONTROL 


RSSF 

VS 

N  (0,1) 


•CONTROL* 


'RSSF* 


NIMEER  OF 


CHI  SQL AR  ED 

TEST 

20  PTS 

30  °TS 

50  PTS 

10? 

3  SC 

421 

498 

5? 

220 

219 

242 

1? 

49 

57 

54 

1C? 

1460 

2407 

3616 

5 % 

1084 

1889 

30  29 

1? 

401 

1011 

1967 

KCLMOG-SMIR 

TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

488 

46  C 

500 

5? 

260 

244 

231 

1? 

64 

50 

28 

10? 

1412 

1782 

2654 

5? 

792 

1083 

1760 

1? 

222 

324 

515 

FOUTZ  FA 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

10? 

503 

494 

5  C 1 

5? 

2  56 

257 

241 

1? 

60 

43 

42 

10? 

2863 

2272 

3962 

5? 

2145 

2678 

3422 

1? 

11  C2 

1575 

2244 

IS  THE  test  OF  UNIF0RM(0,1)  VS  UNIF0RN(C,1) 

IS  THE  RANDOM  STABILIZED  STANDARD  FORM 
ALFFA  =  1.30,  EETA  =  0.0 


HYPOTHESIS  REJECTIONS 


IN  5000  REPLICATIONS 
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CHI  SQUARFO  TFST 


20  PTS 

3C  PTS 

5 C  FT  < 

10? 

3  95 

•ag  c 

5  C  5 

CCATFCL 

5 % 

221 

199 

23  2 

1% 

39 

40 

4  5 

RSS  F 

1C? 

1441 

2116 

3266 

V  S 

53 

10  54 

1598 

2692 

N  (  C  ,  1  ) 

is 

390 

817 

1663 

KOLMOG-  <MIR 

TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

493 

472 

5C2 

CONTROL 

5? 

245 

215 

265 

13 

5C 

47 

52 

PS  S  F 

10? 

1321 

1635 

2362 

VS 

5? 

780 

991 

15C8 

N  (0,1) 

1? 

210 

294 

521 

FOUTZ  FA 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

1C? 

468 

449 

496 

CGNT RCL 

5% 

1  88 

22  2 

237 

1? 

27 

51 

40 

RSSF 

10? 

2479 

2678 

3306 

VS 

53 

18C4 

2059 

27C5 

N  (0,1) 

13 

852 

1080 

1572 

•CON  TROL' 


IS  THE  TEST  0=  UNIFQRM(0,1)  VS  UNIF0RM(C,1) 


*  R  SSF ' 


THE  BANDOM  STABILIZED  STANDARD  FORM 
ALPHA  =  1.60,  EETA  =  0.0 


NLMEER  OF  HYPOTHESIS 


REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  T  EST 
2  C  °TS  30  PTS  50  PTS 


10% 

389 

4C5 

491 

CONTROL 

5  ? 

.  230 

221 

235 

1? 

43 

47 

4  € 

PSSF 

10% 

1326 

1933 

3046 

VS 

5% 

932 

1469 

2428 

N  (0,1) 

1% 

32C 

688 

1440 

KOLMOG-S^IR 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

1C? 

5  27 

474 

489 

CONTROL 

5? 

275 

236 

262 

1% 

59 

38 

54 

RSSF 

10% 

1226 

1521 

2184 

vs 

5? 

760 

9  29 

1389 

t  (0,1) 

1? 

242 

271 

474 

FOUTZ  FN  TEST 

20  PTS  30  PTS  50  PIS 


1  c* 

531 

505 

516 

CCNTFQL 

5? 

265 

263 

2  6  2 

1? 

53 

49 

54 

RSSF 

1C? 

2046 

2274 

2716 

VS 

5? 

1396 

1646 

2081 

N  <C,1) 

1% 

581 

788 

1092 

’CCNTFCL'  IS  1HE  TEST  OF  UNIF0PM(0,1)  VS  UNIF0PM<0,1) 

•RSSF'  IS  TFE  RANDOM  STA8ILIZEC  STANDARD  FOR  N 
ALPHA  *  1.90,  SET  A  =  0.0 


K'PBEF  OF  HYPOTHESIS  REJECTIONS  IN  50  00  R  EP  L  ICA  TI 0N3 
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CONTROL 


RSSF 

VS 

N  (0,1) 


CONTROL 


RSSF 

VS 

N  (0,1) 


CONTROL 


RS  SF 
VS 

N  (0,1) 


'CONTROL' 


•  RSSF' 


M_‘ N B  E  F  OF 


CHI  SQUARED  TEST 

20  PTS  3  C  PTS  50  PTS 


ICS 

4  CO 

457 

496 

5% 

234 

242 

246 

l % 

40 

51 

57 

10  z 

1728 

2935 

4056 

5% 

13  16 

2435 

3629 

1% 

555 

1454 

2612 

KOLMOG— SHIR 

TEST 

20  PTS 

3 C  PTS 

50  PT 

1C* 

504 

474 

477 

5% 

252 

235 

234 

1% 

68 

52 

47 

10* 

1565 

2222 

222  1 

5S 

966 

1410 

2372 

1* 

3  1  C 

481 

901 

FOUTZ  FN  TEST 

20  PTS  3C  PTS  50  FTS 


ICS 

493 

522 

467 

5% 

228 

255 

250 

1? 

49 

53 

46 

10* 

3502 

4014 

4596 

5% 

2859 

2522 

42  9  9 

1% 

1752 

2547 

3487 

S  THE  TEST  0*=  UNIFORMS, 1)  VS  UNIFQRM(0,1) 

IS  THE  RANDOM  STAeiLIZEC  STANDARD  FCRM 
ALPHA  =  l.CO,  EETA  =  0. 25 

HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SCUARED  TEST 


20  PTS  3 C  PTS  50  PTS 


102 

372 

416 

465 

CONTROL 

52 

225 

232 

204 

1% 

42 

56 

45 

RSSF 

10  % 

1521 

2451 

3593 

\S 

5% 

1128 

1933 

3106 

N  (0  ,1  ) 

12 

454 

1057 

2046 

KOLMOG-SMIR 

TEST 

20  FTS 

3  C  PTS 

50  PTS 

102 

486 

512 

483 

CCNTRCL 

52 

242 

272 

2  5  E 

12 

48 

62 

45 

RSSF 

1C? 

1430 

1849 

2  7  Cl 

VS 

5% 

890 

1171 

1882 

N  (0,1) 

U 

268 

271 

6  7  9 

FCUTZ  FN  TEST 

20  PTS  30  PTS  50  PTS 


10? 

461 

465 

455 

CONTROL 

5? 

232 

252 

226 

1? 

37 

53 

47 

RSSF 

10? 

2868 

3284 

40  36 

VS 

52 

22  36 

2684 

3486 

N  (0,1) 

1? 

1174 

1545 

2325 

•CCNTRCL*  IS  THE  T  f  ST  OF  UNI  FORM  (  0 , 1)  VS  UNIFORM  (0,1) 

•RSSF*  IS  TFE  RANDOM  STABILIZED  STANDARD  FORM 
ALPHA  =  1.20,  E  ETA  =  C.25 


N  LM  E  E  F  OF  F  YP  C  7H  E  S  IS  REJECTIONS  IN  5000  REPLICATIONS 
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CFI  SQUARED 

TEST 

20  PT  S 

3  C  PT  S 

50  PT 

1C1 

355 

390 

467 

CONTROL 

51 

223 

212 

2  42 

11 

31 

54 

53 

RSSF 

101 

14  16 

2100 

3277 

VS 

51 

1  C39 

16  10 

27C5 

N  (C,l  ) 

H 

357 

786 

16C7 

KOLMGG-SMIR 

T  EST 

20  RTS 

30  PTS 

5  C  PTS 

101 

514 

518 

5C4 

CCNTFOL 

51 

260 

282 

278 

1  1 

50 

61 

65 

RSSF 

1C1 

1345 

1661 

24C2 

VS 

5? 

807 

1  C17 

1616 

N  (  C  ,1  1 

11 

239 

331 

539 

FOUTZ  FN 

TEST 

2  C  PTS 

3 C  PTS 

50  PTS 

101 

492 

5C2 

495 

CCNTPCL 

51 

248 

260 

24C 

11 

45 

35 

C  - 

RSSF 

101 

23  83 

2650 

3260 

VS 

5? 

16  98 

2023 

26  1  1 

N  10,1) 

11 

822 

1027 

1473 

• CCNTRCL*  IS  THE  TEST  CF  UN!FCRM(0,1)  VS  UNIF0RM<0,1) 

•RSSF*  IS  THE  RANDOM  STABILIZED  STANDARD  FORM 
ALPHA  =  1. 60,  BETA  =  C.25 


MJMeEP  OF  HYPCTHESIS  REJECTIONS 


IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 


20  PTS 

30  PTS 

50  PT 

1C% 

402 

441 

497 

CONTROL 

5  % 

222 

230 

237 

1% 

A3 

45 

46 

RSSF 

VS 

10 % 

1289 

1916 

3  044 

5% 

933 

1458 

2409 

N  (0,1) 

11 

319 

652 

1326 

KCLNCG-SMIR  TEST 


20  PT$  3  C  PTS  50  PIS 


10% 

488 

46  9 

487 

CONTROL 

CO? 

234 

235 

242 

1% 

AO 

52 

C  z 

RSSF 

10% 

1239 

1471 

2111 

VS 

5% 

713 

904 

1357 

N  (0,1) 

1% 

223 

269 

A  A5 

FOUTZ  FN 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

1C3E 

484 

A70 

494 

CONTROL 

5? 

215 

246 

242 

1? 

56 

50 

49 

RSSF 

10% 

2003 

2289 

2773 

vs 

1358 

1688 

2098 

N  <0*1) 

i% 

552 

770 

IOC? 

•CCNTRCL'  IS  THE  TEST  OF  UMFCRM(O.l)  VS  UNIFCRM(0,1) 


•RSSF*  IS  THE  FANDOM  STAEILIZEC  STANDARD  FORM 
ALPHA  =  1,90,  EETA  =  0.25 


MMEEF  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SOU  AR  ED  TEST 

20  PTS  30  PTS  50  FT S 


10? 

3  83 

427 

455 

CC  N  TP  CL 

5% 

213 

248 

253 

1? 

45 

56 

37 

RSSF 

10% 

1732 

2905 

4113 

VS 

5 % 

13  C3 

2395 

3668 

N  (0,1) 

1% 

539 

1432 

2660 

KCLMCG-SMIR 

TEST 

20  PTS 

30  PTS 

50  PTS 

10  * 

5C7 

49  2 

464 

CCNTFCL 

5? 

2  45 

233 

232 

1? 

38 

47 

37 

RSS  F 

1C? 

1663 

2188 

3252 

VS 

5 % 

1011 

1430 

2414 

N  (0,1) 

1? 

289 

455 

93  C 

FOUTZ  FN  TEST 

2 C  PTS  30  PTS  50  PTS 


10? 

507 

492 

534 

CCNTRCL 

5? 

259 

237 

271 

1% 

53 

58 

50 

RSSF 

10? 

34  52 

4011 

4571 

VS 

5% 

2399 

3512 

4297 

N  (0,1) 

It 

1822 

2489 

3465 

•CCNTRCl'  IS  THE  TEST  CF  (JNI  FORM  <0  , 1 )  VS  UNIFCP.M(C,I) 

•RSSF'  IS  THE  FANDOM  STAEILIZED  STANDARD  FORM 
ALPHA  =  I. CO,  0ETA  =  0.50 


NUNS  ER  OF  HYPOTHECS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SOLAR  ED  TEST 


20  PTS  3  C  PTS  50  RTS 


1  c - 

382 

406 

443 

CONTROL 

5% 

213 

19  5 

2  C  7 

1% 

45 

38 

3  6 

RSSF 

10 % 

1665 

2552 

36S5 

V  S 

5% 

1266 

2038 

3222 

N  (0,1) 

li 

509 

i  178 

22C2 

KCLMCG-SMIR 

TEST 

20  PTS 

3  0  PTS 

50  PTS 

10? 

501 

5  03 

486 

CONTROL 

5? 

248 

256 

2  53 

1? 

45 

55 

57 

PSSF 

10? 

16CC 

2092 

30  2  5 

VS 

5? 

1026 

1425 

2208 

N  (0,1) 

1? 

342 

532 

951 

FOUTZ  FN  TEST 

20  PT S  3 C  PTS  50  °TS 


1C? 

5  06 

467 

471 

CONTROL 

5? 

246 

236 

228 

1? 

44 

46 

43 

RSSF 

10? 

28  16 

2278 

3920 

VS 

5? 

2168 

2680 

3382 

N  (0,1  ) 

1? 

1180 

1546 

2270 

•CCMPCL*  IS  11-  E  TEST  OF  UMFCRM<0,1)  VS  UNIFORMS, 1) 

'  RSSF*  IS  THE  PANDCM  STAeiLIZEC  STANDARD  FORM 
ALPHA  =  1.30,  EETA  =  0. 50 


M.PBEF  CF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 


CONTROL 


RSSF 

VS 

N  (0,1) 


CONTROL 


PSSF 

VS 

N  (0,1) 


CONTROL 


RSSf 

VS 

N  (0,1) 


•CONTROL' 


•  RSS  F* 


NLMEER  OF 


20  PTS 

3  C  PTS 

50  PTS 

1C? 

413 

440 

470 

5  ? 

251 

222 

2C5 

1% 

36 

48 

41 

1  C? 

1451 

2269 

3435 

5? 

1063 

1794 

2841 

1? 

366 

948 

1626 

KCLNCG-SMI  R 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

10? 

478 

518 

513 

5? 

2  52 

278 

277 

1? 

49 

65 

6C 

10? 

1435 

1853 

2636 

5? 

888 

1229 

1881 

1? 

2  67 

444 

758 

FOUTZ  FN 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

1C? 

493 

463 

488 

5? 

241 

220 

227 

1? 

47 

47 

40 

10? 

2314 

2643 

3268 

5? 

16  86 

202  8 

2594 

1? 

761 

1 C52 

1449 

S  THE  TEST  OF  UMFCRM(O.l)  VS  J  NI  FORM  (0,  1  ) 


IS  THE  RANDOM  ST  A  El  LI 
ALPHA  =  1.60,  EETA 


HYPOTFES  IS  REJECTIONS 


EC  STANDARD  FORM 
=  0*50 


IN  5000  R  EP  LI  CAT  I  C  NS 
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CHI  5QU  AR  EO  TEST 


20  PTS  30  PTS  5 C  PTS 


1C? 

419 

411 

466 

CCNTFOL 

5 % 

218 

230 

234 

1? 

38 

54 

52 

RSSF 

1C? 

1299 

1934 

3  C42 

VS 

5? 

9  16 

1438 

24C7 

N  (  C,1  ) 

1% 

297 

680 

1348 

KOLMOG-  SMIR 

TEST 

20  PTS 

30  PTS 

50  PTS 

10% 

5  C2 

494 

48  2 

CONTPCL 

5% 

26C 

255 

248 

1? 

47 

58 

62 

RSSf 

1C? 

1239 

1559 

2146 

VS 

5? 

738 

961 

1357 

N  (C,l) 

1? 

214 

296 

475 

FCUTZ  FN 

TEST 

2  0  PTS 

30  FTS 

50  PTS 

10% 

473 

479 

508 

CONTROL 

5% 

259 

229 

240 

1% 

5  5 

53 

52 

RSSF 

10? 

2018 

2310 

2708 

VS 

5? 

1384 

1629 

2043 

N  (0,1) 

1? 

591 

749 

1023 

•CONTROL'  IS  THE  TEST  OF  UNIFORM(0,1>  VS  UNIFORMOtl) 

*  R  S  SF*  IS  THE  FANDOM  STAEILIZEO  STANDARD  FORM 
ALPHA  =  i.90,  E  ETA  =  0.50 


MM8EF  OF  FYPCTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 


20  PTS  3 C  PTS  50  PTS 


1C? 

410 

445 

436 

CONTROL 

c? 

229 

236 

221 

1? 

30 

49 

46 

RSSF 

10? 

1831 

2913 

4 1  C  8 

VS 

5  % 

13  58 

239  1 

3655 

b  (C,li 

1% 

566 

*  43  8 

2690 

KOLMOG-SMIR 

T  EST 

2  C  FT  S 

30  PTS 

5  C  PTS 

10? 

524 

505 

527 

CCNTfQL 

5? 

250 

259 

25  I 

1? 

46 

53 

55 

PSSf 

1C? 

1627 

2196 

33C4 

V  s 

5% 

10  51 

1421 

2387 

N  10,1) 

1? 

3  04 

500 

931 

F CUT Z  FN  TEST 


2  C  PTS  30  PTS  50  PTS 


10? 

515 

497 

502 

CONTPCL 

5? 

229 

254 

2  c  c 

1% 

51 

59 

46 

PSSF 

10 % 

3586 

4057 

4542 

VS 

5% 

3005 

2591 

4256 

N  <0,i) 

1? 

1893 

2495 

3453 

*  CCNTRCL*  IS  THE  T£$T  CF  UNIFORMS, 1)  VS  UNIFORM<0,1) 

'  R  S  SF  *  IS  THE  RANDOM  STABILIZED  STANDARD  FORM 
ALPH;  =  1.00,  EETA  =  C.75 


NUMBER  OF  HTPCTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SOUAREO  TEST 


20  PTS  2C  PTS  50  PTS 


1C? 

368 

447 

505 

CCNTPOL 

5% 

2  CC 

241 

222 

1? 

43 

58 

39 

PSSF 

101 

19  41 

2511 

4  C97 

VS 

51 

1495 

2429 

36  T  C 

N  ( C,l) 

1? 

704 

1535 

2762 

KQLMOG-SMIR 

TEST 

20  PTS 

3C  PTS 

50  PTS 

10? 

4  75 

50  4 

524 

CCNTPCL 

5? 

242 

237 

264 

1? 

47 

45 

47 

RSSP 

10? 

20  17 

2  64C 

3726 

VS 

51 

1331 

1930 

3089 

N  <C,1) 

1? 

546 

904 

1762 

FCUTZ  FN  TEST 

20  PTS  30  PTS  50  PTS 


10? 

461 

513 

490 

CCNTPOL 

5  1 

226 

248 

248 

1 % 

45 

54 

5  5 

RSSF 

1C? 

2895 

2212 

4  C  5  1 

VS 

5? 

2227 

2712 

3542 

N  (0  ,1  ) 

11 

12C9 

1595 

2426 

'CCNTPCL  *  IS  THE  TEST  CF  UNIFORM(Otl)  VS  UNIFCPMC,1) 

•PSSF*  IS  THE  RANDOM  STAEILIZEC  STANDARD  FORM 
ALPHA  =  1.20,  EET  A  =  C.75 


HUMBER  CF  HYPOTHESIS  REJECTIONS  IN  5GC0  REPLICATIONS 
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CHI  SQUARED  TEST 

20  PTS  3  C  PTS  50  PTS 


1C? 

4  CO 

407 

501 

CONTROL 

5? 

2  43 

225 

256 

1  ? 

48 

51 

46 

RS  SF 

10? 

1628 

2451 

3633 

VS 

5? 

1240 

15  6  5 

3C99 

N  (0,1) 

1? 

509 

1126 

21 1  C 

KOLMOG-SNIR 

TOST 

20  PTS 

3 C  PTS 

50  FTS 

1C? 

491 

505 

476 

CCNT  FOL 

111 

244 

243 

1? 

48 

54 

50 

F  S  SF 

10? 

1622 

2  IC4 

302  3 

VS 

5? 

10  50 

1465 

2330 

b  (0,1) 

1? 

375 

588 

1145 

FOUTZ  FN  TEST 

20  PTS  30  PTS  50  FTS 


1C? 

522 

490 

45  C 

CONTROL 

5? 

266 

26  7 

243 

1? 

52 

60 

56 

RS  SF 

10? 

2426 

2722 

3367 

VS 

5? 

1753 

20  67 

2727 

N  (0,1) 

1? 

8  25 

1C58 

1614 

•CONTROL*  IS  THE  TEST  OF  UNIFORM(Otl)  VS  UNIFORM(Otl) 

•RSSF*  IS  THE  RANDOM  STAEILIZEC  STANDARD  FORM 
ALPHA  =  1 . 60 1  3ETA  =  0.75 


M'NBEF  OF  HYPOTHESIS  REJECTIONS  IN  5000  R?d  LI  CAT  I  C  HS 


CHI  SCU6RED  TEST 


2  C  PTS 

30  PTS 

50  PT 

10* 

405 

4C7 

469 

CONTROL 

5* 

?38 

229 

215 

1* 

54 

49 

46 

FSSF 

10* 

1366 

1989 

3088 

VS 

5* 

987 

1450 

2487 

N  10,1) 

1* 

341 

704 

141  S 

KOLVOG- 

SMIR 

TEST 

20  PTS 

3  C  PTS 

50  PT 

CONTROL 

10* 

50  1 

505 

483 

5* 

246 

268 

245 

1* 

57 

53 

46 

PSSF 

10* 

1273 

1558 

2246 

VS 

5* 

797 

936 

1465 

N  (0,1) 

1* 

234 

308 

508 

FOUTZ  F 

N  TEST 

20  PTS 

3C  PTS 

50  FT 

1C* 

497 

457 

478 

CONTROL 

5* 

242 

227 

242 

1* 

58 

48 

48 

PSSF 

10* 

2102 

2256 

2745 

VS 

5* 

1479 

16  15 

2  C  E 1 

N  (C,l) 

1* 

647 

76  0 

1076 

i 

;CNTROL*  IS  THE 

TEST 

OF  UNI 

FCRM< 

:o,i)  VS  U 

NI FORM  (0 

•  RSS  F' 

IS  THE 

RANDOM  STABILIZED  STANDARD  FCRM 

ALPHA  = 

1 »  SC, 

EETA 

=  0.75 

'LM9EF  OF 

HYPCTH 

ESIS 

REJECT 

IONS 

IN  5000  R 

EPLICATI 
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CHI  SQUARED  TEST 


20  PTS  3  C  PTS  50  PTS 


1C? 

3  55 

441 

455 

CONTROL 

5% 

213 

238 

2  43 

1% 

35 

53 

62 

PSSF 

10? 

1793 

2  52  C 

4109 

VS 

5? 

1355 

2  38  6 

3658 

1*  (Otl) 

1? 

554 

1413 

2707 

KGLMOG-SMIR 

TEST 

20  PTS 

20  P^S 

50  FT 

10? 

493 

504 

457 

CONTROL 

5% 

241 

231 

228 

1% 

50 

42 

43 

RSSF 

1C? 

1640 

2203 

3325 

VS 

5% 

10  04 

139  2 

2434 

N  (  C  ,1  ) 

1% 

295 

421 

935 

FOUTZ  FN 

TEST 

20  PTS 

30  PTS 

50  PT 

10? 

526 

488 

48  1 

CCNTRCL 

5? 

260 

247 

223 

I? 

54 

58 

50 

RSS  f 

1C? 

3530 

4001 

4577 

VS 

5? 

29  33 

3501 

4284 

N  ( C  »1) 

1? 

1849 

2405 

34C5 

•CONTROL'  IS  THE  TEST  OF  UNIFORM(0,1)  VS  UNI  FORM  (0 , i ) 

•  FSSF  •  IS  THE  RANDOM  STABILIZED  STANCARD  FORM 
ALPHA  =  1 .00  t  BETA  =  1.00 


NOME  EP  QF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  S  QU  AR  £0  TEST 

20  PTS  2C  PTS  50  PTS 


1  Cl 

363 

428 

518 

CONTROL 

5? 

218 

226 

242 

II 

45 

44 

55 

RSSF 

1C? 

2651 

2826 

4682 

VS 

5? 

2291 

3457 

4472 

N  (  C  ,1  ) 

1? 

1344 

2572 

3935 

KOLMOG-SMIR 

T  EST 

20  PTS 

30  PTS 

50  PTS 

10? 

526 

533 

5  24 

CONTROL 

5? 

2  43 

263 

2  62 

1? 

51 

48 

7C 

PSS  f 

10? 

2937 

2776 

4636 

vs 

5? 

23  26 

3235 

4378 

N  <C,1) 

1? 

1245 

2035 

3452 

F  CUT Z  F N  TEST 

20  PTS  30  PTS  50  PTS 


10? 

488 

477 

500 

CONTROL 

5? 

259 

216 

236 

1? 

54 

53 

44 

PSS  F 

10? 

2962 

3429 

4162 

VS 

5? 

2260 

2813 

3660 

N  tO,l) 

1? 

1237 

170  5 

25  50 

'CONTROL  *  IS  THE  TEST  OF  UNIF0RM(0,1)  VS  UNIFORMS, 1) 

*  RSSF'  IS  TFE  FANDOM  STAEILIZEC  STANDARD  FORM 
ALPHA  =  1.  20,  E  ETA  =  1  *00 


MMEEP  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED 

TEST 

20  PTS 

3C  PTS 

50  FT 

101 

412 

443 

472 

CCNT FOL 

51 

220 

237 

220 

11 

30 

50 

46 

RSS  f 

1  Cl 

1868 

2768 

3  5  5  8 

V  S 

51 

1487 

2294 

3541 

N  (  C  ,1  ) 

11 

701 

1427 

2578 

KCLMOG-SMIR 

TEST 

20  PTS 

3C  PTS 

50  PT 

101 

513 

477 

5  14 

CONTROL 

51 

268 

251 

251 

11 

56 

56 

49 

RSSF 

101 

1968 

2603 

3652 

VS 

51 

1370 

1903 

30  16 

N  (0,1) 

11 

586 

905 

1765 

FQUTZ  F  fi  TEST 

20  PTS 

3 C  PTS 

50  PT 

101 

5  07 

504 

5C5 

CONTROL 

51 

25C 

240 

252 

11 

58 

45 

58 

RSSF 

101 

2414 

2792 

3409 

VS 

51 

1786 

2120 

2771 

N  (C,l) 

11 

853 

1123 

1625 

*  CONTROL '  IS  THE  TEST  CF  UNIF0RM<0,1)  VS  UNIFORMS, 1) 


*  RSSF*  IS  THE  RAN COM  STABILIZED  STANDARD  FO«M 
ALPHA  =  1.60,  EETA  =  1  .OC 


NI.MEER  OF  HYP  CTH  E  SIS  REJECTIONS 


IN  5000  REPLICATIONS 
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CHI  SQLARED 

TEST 

20  FTS 

3  C  PTS 

50  PT 

ic* 

3  74 

431 

525 

CONTROL 

5? 

231 

232 

2  6  3 

1% 

48 

57 

58 

FSSF 

10% 

1329 

15  53 

3  1C  7 

VS 

5% 

977 

1438 

2466 

N  (0,1) 

1% 

347 

708 

1459 

KQLMOG-SMIR 

TEST 

20  PTS 

3 C  PTS 

50  FTS 

1C? 

494 

48  1 

516 

CCNTPOL 

5 X 

?-3g 

240 

241 

1% 

"54 

47 

43 

RSSF 

10? 

1322 

1625 

23  C  c 

VS 

5  % 

8C8 

978 

1533 

N  (C,U 

1% 

225 

3  02 

540 

FOUTZ  FN 

T  EST 

20  FTS 

30  PTS 

5 C  PTS 

10? 

525 

523 

520 

CCNTFCL 

5? 

236 

271 

296 

1? 

50 

67 

55 

RSSF 

1C? 

2086 

2  280 

2758 

VS 

5% 

1409 

1645 

2075 

N  (  C  ,1  ) 

1? 

620 

71  5 

1078 

*C  CNT  F  CL  '  IS  THE  TEST  OF  UNIFQRM(C,1)  VS  UMPCPM(0,1) 

'FSSF  '  IS  THE  FANCCM  STABILIZEC  STANDARD  FORH 
ALPHA  =  1  .90  t  EETA  =  1 .00 


N  L  H  6  E  P  CF  HYPOTHESIS  REJECTIONS  IN  50CO  REPLICATIONS 
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CHI  SQUARED  TEST 


CCNTPGL 

N  IXNCRV 

vs 

N  (  C ,1  ) 


CONTROL 


MIXNCRM 

VS 

N  (C,l) 


CCNTROL 


MXNCPM 

VS 

N  (0,11 


2 C  PTS 

3 G  PTS 

50  PTS 

10? 

388 

469 

4  £6 

5? 

2  29 

238 

224 

1? 

39 

48 

43 

1C? 

1272 

1925 

30C4 

5? 

9  09 

1407 

2381 

1? 

3  29 

664 

1337 

KOLMOG-SMIR 

TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

52  5 

513 

5  10 

5? 

2  54 

260 

2  4  S 

1? 

42 

46 

64 

10? 

1235 

1537 

2133 

5? 

7  13 

S3 1 

1369 

1? 

208 

28  5 

44  S 

FOUTZ  FN 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

10? 

508 

476 

523 

5? 

245 

222 

255 

1? 

50 

38 

6C 

10? 

1982 

2191 

26  10 

5? 

1336 

1576 

1955 

1? 

533 

672 

977 

•CCNTROL*  IS  THE  TEST  CF  UNIFQRM(0,1)  VS  UNIFORMOtl) 

•MXNCRM  IS  TEE  MIXED  NORMAL  DISTRIBUTION 
C.  C  *N  (.0, 1)  +  1  .OC*N<  C,2  ) 


rUNBEP  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATION 
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CHI  SQUARED  TEST 


20  PTS  30  PTS  5  C  FTS 


10? 

4CS 

411 

458 

CCNTFCL 

5? 

233 

205 

234 

1? 

29 

34 

46 

f'lXNCRM 

10? 

2513 

3558 

462  C 

VS 

5* 

2088 

3C8  0 

4342 

N  (0,1) 

12 

1028 

2052 

36  C  5 

KOL  NOG-  SMI  R 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

10? 

475 

463 

50C 

CONTROL 

5? 

238 

229 

242 

1? 

46 

35 

38 

NIXNCPM 

10? 

21  88 

2827 

39  6  5 

VS 

5  ? 

1420 

1974 

3161 

N  (0*1) 

1? 

512 

727 

1515 

FOUTZ  F  1^  TEST 

20  PTS  3 C  PTS  50  PTS 


10? 

522 

454 

481 

CONTROL 

5? 

236 

217 

226 

1? 

47 

41 

52 

MXNCFM 

10? 

3338 

3712 

4346 

VS 

5? 

2724 

3174 

3935 

N  (0,1) 

1? 

1657 

2064 

2970 

'  CONTROL'  IS  THE  TEST  CF  UNIFORMS, 1)  VS  UN!F0RM(C,1) 

•MXNCRM*  IS  THE  MIXEC  NORMAL  DISTRIBUTICN 
C.C  *M.0,1)  ♦  1.00*N<0,3) 


MJNBEF  OF  HYPOTHESIS  REJECTIONS  IN  5000  R  EP  L  IC  A  TIONS 
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CHI  SQUARED  TEST 

20PTS  3CPTS  50PTS 


10? 

407 

425 

5  11 

CONTROL 

5% 

2  34 

232 

252 

1% 

48 

4  £ 

5C 

MIXNOFM 

10? 

3427 

4423 

4925 

VS 

5% 

3006 

4154 

4863 

N  (0  ,1) 

1? 

1  7  £1 

3377 

46  29 

KOLMOG-SMIR 

TEST 

20  PTS 

3C  PTS 

50  FTS 

1C? 

510 

490 

511 

CCNTFOL 

5% 

261 

241 

250 

1 % 

59 

48 

47 

MXNCRM 

1C? 

2976 

2799 

4717 

VS 

5? 

2113 

3025 

42  £6 

N  (  Ctl  ) 

1? 

890 

1471 

2792 

FOUTZ  FN  TEST 

20  FT S  30  PTS  5C  PTS 


10? 

518 

501 

528 

CONTROL 

5 % 

242 

268 

276 

1% 

68 

66 

62 

MIXNCRM 

1C? 

4131 

4502 

4  £6  5 

VS 

5? 

2686 

4195 

4742 

N  (C,l) 

1? 

2653 

2330 

4272 

•CCNTRCC*  IS  THE  TEST  OF  UN I^ORM ( C  ,1 )  VS  UNIFCRM<C,1> 

•MIXNCRM'  IS  THE  MIXED  NORMAL  DISTRIBUTION 
0.0  *N<  .0»1  )  +  1  »00*N  ( 0  »4  1 


NUMBER  OF  HTPCTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SOU  £R ED 

TEST 

20  PTS 

20  PTS 

5  C  PTS 

10? 

3  82 

422 

510 

CCNT POL 

5 % 

215 

137 

254 

1% 

41 

38 

46 

HmCRH 

10? 

472 

542 

713 

vs 

5? 

2  53 

308 

3  70 

F  (Ctl) 

1% 

29 

66 

74 

KOLMOG-SMI  R 

TEST 

20  PTS 

30  PTS 

50  PTS 

1C? 

488 

537 

525 

CCNTFCL 

5% 

238 

250 

265 

1? 

46 

50 

53 

M  IXNCRM 

1C? 

587 

667 

655 

VS 

5? 

287 

328 

367 

h  ( c  a  ) 

1? 

54 

70 

84 

FOUTZ  FN 

TEST 

20  FT  S 

30  PTS 

50  PTS 

10? 

468 

480 

517 

CCIVTPCI 

5? 

220 

258 

247 

1? 

46 

64 

56 

MIXNCRN 

1  C? 

799 

796 

8  El 

VS 

5? 

425 

426 

477 

N  (oai 

1? 

105 

107 

136 

•  C  CNTPC L '  IS  THE  TEST  CF  UNIFORMS, 1)  VS  UNIF0RM(C,1) 

'  M IXNQRM  *  IS  THE  MIXED  IVCPMAL  DISTRIEUTION 
C.7C*N(.0,1>  +  C.30*N(0,2> 


M.MEEP  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATICNS 
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CHI  SQUARED  TEST 


20  PTS  3C  PTS  5 C  PTS 


10  ? 

421 

457 

510 

CCNT FOL 

5? 

245 

244 

22  1 

1% 

49 

52 

3  € 

MIXNCRM 

1C? 

551 

727 

1C4  7 

VS 

5? 

345 

438 

602 

N  (  C  ,1  ) 

1% 

82 

131 

185 

KOLMCG-SMIR 

TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

490 

504 

527 

CCNTROL 

5? 

251 

256 

23  S 

1? 

46 

58 

42 

MIXNORM 

10? 

6  46 

745 

845 

VS 

5? 

337 

405 

447 

N  (0.1) 

1? 

81 

93 

E3 

FQUTZ  F I'  TEST 

20  PTS  3  C  PTS  50  PTS 


10? 

5  19 

494 

507 

CONTROL 

5? 

24C 

242 

256 

1? 

55 

46 

51 

MIXNCRM 

10? 

1029 

1139 

1249 

VS 

5? 

59  1 

671 

778 

h  (0,1) 

1? 

172 

193 

262 

•CCNTRCL*  IS  TFE  TEST  CF  UNIFORMS, 1)  VS  UNI  FOR  V  (  Ct  I ) 

•MIXNORM  IS  THE  MIXED  NORMAL  DISTRIBUTION 
0.7  C*N  (.0,1)  +  0.30*N(0»3) 

KIMBEF  CF  HYPCTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SOU  JR  EO  TEST 


20  PTS  3  C  P7S  50  PTS 


1C? 

400 

435 

5  15 

CONTROL 

5? 

237 

223 

240 

1 ? 

45 

50 

52 

M  I  >  N  CR  M 

10? 

6  36 

910 

1320 

VS 

1% 

4C5 

564 

336 

N  (0,1  ) 

1% 

109 

175 

222 

KOLMOG-SMIR 

TEST 

2 C  PTS 

30  PT$ 

50  PT 

10? 

508 

531 

5  C  E 

CCNTFCL 

5% 

259 

259 

248 

1% 

54 

51 

54 

N  IXNCRM 

1C? 

747 

853 

10C1 

VS 

5? 

407 

441 

532 

N  (C,l  ) 

1? 

1  15 

122 

148 

FO'JTZ  FN  TEST 

20  FT 5  30  PTS  50  PTS 


10? 

4  92 

478 

511 

CONTROL 

5  % 

2  29 

232 

241 

1? 

61 

47 

42 

MIXNCRM 

10? 

1281 

1418 

16  66 

vs 

5? 

774 

888 

1121 

N  (0,1) 

1? 

238 

222 

4  C 1 

*  CCNTRC  L*  IS  THE  TEST  OF  UNI  FOR  M  (  0  ♦  1 )  VS  UNIFORMS, 1) 

' M  I XNORM •  IS  THE  MIXED  NCRMAL  DISTRIBUTION 
C.7C*N(.0,1)  +  0.30*N<0,4) 


MMEEP  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 


73 


CHI  SQU £R  ED  TEST 


20  PTS 

30  PTS 

50  FT  S 

10? 

3  80 

431 

5C7 

CCNTFCL 

53 

221 

245 

220 

13 

40 

59 

56 

MXNCRM 

1C? 

388 

484 

545 

V  S 

5? 

214 

274 

275 

N  <G,1) 

1? 

40 

63 

62 

KOLMOG- 

-  SM  IR  TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

468 

480 

48  1 

CCNTRCL 

5? 

220 

236 

232 

1? 

53 

46 

28 

MIXNCRM 

10? 

519 

539 

556 

VS 

5? 

258 

277 

286 

N  (0,1) 

1? 

58 

55 

53 

FOUTZ 

F  It  TEST 

20  PTS 

3  C  °  T  S 

50  PTS 

10? 

487 

498 

502 

CONTROL 

5? 

234 

242 

265 

1? 

42 

45 

54 

VIXNC  FM 

10? 

684 

689 

687 

VS 

5% 

352 

354 

400 

N  (0,1) 

1? 

83 

92 

96 

•CONTROL*  IS  THE  TEST  CF  UNIP0RM(0,1)  VS  UNIFORM(C,l) 

•NIXNCPM*  IS  THE  MIXED  NORMAL  DISTRIBUTION 
C.  €C*N  (  .0*  I )  +  0.20*N(0,2) 


t'UMBEF  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATICNS 


74 


CHI  SQUARED  TEST 

20  PTS  30  PTS  5 C  PTS 


1C? 

417 

405 

*92 

CONTROL 

5? 

218 

211 

242 

1? 

43 

38 

42 

b  IXNCRM 

1C? 

477 

544 

736 

VS 

5? 

281 

306 

364 

N  ( C,1  ) 

1 % 

48 

64 

58 

KOLMOG-SMIR 

T  EST 

20  FTS 

30  PTS 

50  PT 

CCNTRCL 

10? 

5  CC 

469 

527 

5? 

267 

222 

2  6  C 

1? 

47 

46 

48 

MIXNCRM 

1C? 

587 

590 

678 

vs 

5? 

3  17 

297 

342 

N  *  (  C  »1 ) 

1? 

69 

59 

78 

FCUTZ  FN  TEST 

20  PTS  30  PTS  50  PTS 


10? 

495 

46? 

479 

CONTROL 

5? 

234 

213 

245 

1? 

49 

44 

48 

NIX  NC  PM 

10? 

839 

819 

945 

VS 

5? 

475 

454 

522 

N  (Ctl) 

1? 

1  18 

130 

142 

’CCNTRCL '  IS  THE  TEST  OF  UNIF0RM(0»1)  VS  UN  IFQR  N  <  0 , 1 ) 

•MIXNCRN'  IS  THE  MIXED  NORMAL  DISTRIBUTION 
C  •  £C*N  <•  Ot  1)  +  0.2C*N<0,2) 


NUMBER  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 


20  PTS 

3  C  PTS 

50  PTS 

1C* 

402 

4  50 

539 

CCMFCL 

5? 

225 

242 

247 

1* 

43 

49 

6C 

MIXNGRM 

1C* 

5  07 

637 

9C2 

\S 

51 

310 

361 

496 

N  <C,1) 

11 

74 

89 

138 

KOLMOG-  SMIR 

TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

511 

502 

542 

CONTROL 

51 

253 

252 

270 

1? 

51 

43 

63 

NIXNCRM 

10? 

622 

661 

712 

VS 

5? 

3  29 

360 

416 

N  (0,1) 

1? 

85 

69 

93 

FCUTZ  FN 

TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

485 

496 

514 

CONTROL 

5? 

238 

243 

259 

1? 

48 

48 

c  c 
<0*  +* 

M  IXNQRM 

10? 

958 

1011 

1168 

VS 

5% 

548 

583 

721 

N  (0,1) 

1? 

171 

174 

246 

•CONTROL*  IS  THE 


TEST  CF  UNI  FCRM  (0*1)  VS  UNIFORMO,!) 


•  M  1  >NCR  V’  IS  TFE  MIXED 
C.£C*Nl.O,l)  + 


NORMAL  D  1ST  RI8UTI  ON 
0.2  0*N<  C,4  ) 


H  NB  E  F  CF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATION 


76 


CHI  SQUARED  TEST 

20  PTS  3 C  FT S  50  PTS 


1C3 

420 

470 

47C 

control 

53 

231 

240 

236 

1 3 

40 

51 

44 

MI >NCPM 

10? 

4 1  C 

478 

525 

VS 

53 

228 

241 

247 

N  (0,1) 

13 

35 

48 

46 

KOLMOG-SNIR 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

CONTROL 

1C? 

499 

490 

519 

53 

244 

245 

251 

13 

54 

44 

55 

MXNCRM 

10? 

521 

512 

535 

VS 

5? 

264 

257 

266 

N  <C,1> 

1  3 

57 

47 

64 

FOUTZ  FN  TEST 

20  PTS  30  PTS  50  PTS 


1C? 

493 

490 

4  £  C 

CONTROL 

53 

245 

256 

229 

13 

43 

62 

47 

HXNCRN 

10? 

577 

,  583 

587 

\S 

53 

282 

310 

2  S  C 

N  (C,l) 

13 

66 

73 

77 

•CONTROL*  IS  THE  TEST  OF  UNIFCRMO.l)  VS  JNIFQRMOtl) 

•MIXNORM*  IS  THE  MIXED  NORMAL  E ISTP.  IBUT  I CN 
0.907NI  .0,1.)  +  0.10*N(0,2) 


NUME  EP  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 


77 


CHI  SQU  £R  ED  TEST 

20  PTS  30  PTS  50  FT S 


10? 

432 

416 

4  £6 

CONTROL 

5? 

246 

214 

215 

1? 

43 

41 

4  C 

V  IXNCRM 

10? 

4  56 

431 

537 

VS 

53! 

262 

233 

2  58 

n  cc,n 

1? 

51 

43 

53 

KOLMOG-SMI  R 

T  EST 

20  PTS 

30  PTS 

50  PTS 

10* 

5  C9 

495 

509 

CCNTFOL 

5? 

271 

239 

25  2 

1? 

55 

46 

55 

MIXNCRM 

1C? 

553 

5-26 

557 

VS 

5? 

302 

26  1 

274 

N  (C  ,1  ) 

1? 

67 

53 

64 

FOUTZ  FN  TEST 

20  PTS  3 C  PTS  50  PTS 


10? 

516 

471 

520 

CCNTFOL 

5? 

271 

239 

252 

1? 

58 

51 

47 

MIXNCRM 

10? 

66  C 

653 

675 

VS 

5? 

353 

357 

356 

N  (0,1) 

1? 

95 

88 

£  5 

•CONTROL'  IS  THE  TEST  CF  UNIF0RM<0, 1)  VS  UNIFCRM(0,1J 

•NI>NCPy  IS  THE  MIX  EC  NORMAL  DISTRIBUTION 
C*5C*N(.0,1)  +  0.1C*N(0,3) 


NUMEER  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 

20  PTS  30  PTS  50  PTS 


1C? 

357 

439 

505 

CON  TROL 

5? 

240 

244 

246 

1? 

31 

53 

55 

R  IXNCRN 

1C? 

419 

50  8 

551 

VS 

5? 

260 

266 

3C2 

N  (C,l) 

1? 

36 

65 

74 

KCLMQG-SMIR 

T  EST 

2  C  PTS 

3C  PTS 

5  C  PT 

CCNTfOL 

10? 

4  84 

516 

48  £ 

5? 

240 

247 

229 

1? 

56 

47 

45 

HIXNCRN 

1C? 

525 

534 

54  £ 

V  s 

5? 

268 

28  2 

27  7 

N  (0,11 

1? 

64 

61 

6  C 

FO'JTZ  FN  TEST 


2  C  PTS  30  PTS  50  PTS 


10? 

514 

511 

45  1 

CCNTRCL 

5? 

7  47 

253 

245 

1? 

48 

50 

58 

NIXNCRM 

10? 

736 

707 

747 

VS 

5? 

385 

384 

422 

N  (0,1) 

1? 

81 

105 

112 

•CCNTRCL'  IS  THE  TEST  CF  UNIFORMS, 1)  VS  UNIFCRM(C,1) 

•RIXNCRN'  IS  THE  HIXEC  NCRMAL  CISTRIBUTION 
C. 90*N (  .0,1  )  +  0.10*N<0,4) 


NUMBER  CF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED 

TEST 

20  PTS 

30  PTS 

50  PTS 

10% 

3  €6 

417 

494 

control 

5% 

225 

241 

238 

1% 

39 

44 

46 

R  IXNCRR 

10% 

1459 

2008 

3146 

VS 

5% 

1  093 

15  3  3 

2556 

N  (C,l) 

1% 

407 

793 

1568 

KCLMOG-SMIR 

TEST 

20  PTS 

3C  PTS 

5  C  PTS 

10% 

493 

483 

487 

CCNTFOL 

5% 

267 

246 

238 

1% 

52 

53 

36 

MIXNCRR 

1C% 

2179 

2806 

3756 

VS 

5% 

16C6 

2181 

3226 

N  <C,1) 

1% 

714 

1148 

2  043 

FOUTZ  FN 

TEST 

2  C  PTS 

30  PTS 

50  PTS 

10% 

535 

495 

499 

CONTROL 

5% 

250 

246 

248 

1% 

54 

43 

41 

MIXNGRM 

10% 

18  09 

1994 

2492 

VS 

5% 

1198 

1428 

1825 

N  (0,1) 

1% 

509 

60  2 

865 

’CCNTRCL'  IS  THE  TEST  CF  UNIFORMED  ,11  VS  UNIFORM  (0,1) 

*  M  I  >NO  PM '  IS  THE  MIXEO  NORMAL  DISTRIBUTION 
C»  7 C*N ( .5, 1)  +  0  .3C*N( C,3) 


NUMBER  OF  HYPOTHESIS  REJECTIONS  IN  5000  REP  LICATI CHS 
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CHI  SQUARED  TEST 

20  PTS  30  PTS  50  FT S 


102 

377 

449 

512 

CCNTPCL 

5? 

20  5 

239 

257 

IS 

46 

46 

44 

NIXNCRM 

1C? 

1557 

2112 

33C4 

VS 

5? 

1170 

1563 

2660 

N  10,1) 

1? 

473 

750 

1 6  4  C 

KOLf'CG-SMI  R 

TEST 

20  PTS 

30  PTS 

50  PTS 

10? 

482 

481 

538 

CONTROL 

5? 

229 

260 

271 

IS 

45 

69 

56 

MIXNCRM 

10? 

24  5  C 

3101 

4056 

VS 

5? 

1855 

2510 

3596 

N  (0,1) 

1? 

9  Cl 

1344 

24CQ 

FOUTZ  FN  TEST 

20  PTS  3  C  PTS  50  PTS 


1C? 

479 

512 

531 

CONTROL 

5? 

233 

262 

255 

1? 

50 

53 

53 

MXNCRN 

10? 

1781 

1917 

2405 

VS 

5? 

12C9 

1369 

1774 

N  (0,1) 

1? 

462 

573 

842 

•CCNTRCL'  IS  THE  TEST  CF  UNI  FORM (0 » 1)  VS  UN  I FCR M ( C , I  ) 

•  ^  I X  N  CR  N  •  IS  THE  MIXEC  NORMAL  CISTP.IBUTICN 
0.£C*M.5,1>  +  0 . 20*  N(  0,3) 


NUMBER  CF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 
2 C  PTS  3C  PTS  50  PTS 


10* 

389 

439 

490 

CCNTFOL 

5* 

221 

251 

221 

1* 

38 

55 

45 

f- IX  N  CRM 

1C? 

1568 

2228 

3377 

VS 

5? 

1177 

1736 

2832 

N  (  C,1  ) 

1? 

507 

551 

1756 

K0LMOG-  SMR 

TEST 

20  PTS 

30  PTS 

50  PT 

10* 

535 

489 

490 

CCNTRCL 

5? 

258 

253 

245 

1* 

55 

54 

35 

VIXNCRM 

10* 

2467 

2194 

41C5 

VS 

5? 

1865 

2553 

3659 

N  (C,l) 

1% 

857 

1438 

2511 

FCUTZ  FN  TEST 

20  PTS  30  PTS  50  PTS 


10? 

502 

456 

495 

CCNTRGL 

5? 

2  52 

250 

235 

1? 

62 

51 

MXNCFM 

10? 

183  8 

2156 

2665 

VS 

5? 

1220 

1551 

1978 

N  (0,1) 

1? 

486 

701 

975 

’CONTROL'  IS  THE  TEST  CF 


UNI  FORM ( 0  »1 )  VS  UN  IFOR  MO  1 1 ) 


•KI>NCRf»<  IS  TFE  MIXED 
C.  EC*N  («5»  1)  + 


NORMAL  DISTRIBUTION 
0.2  0*N(  C,4  ) 


NL'FBEF  CF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 


CHI  SQUARED 

TEST 

20  PTS 

3  C  D  T  S 

50  PTS 

10? 

431 

413 

5  2  6 

CCNTFCL 

5? 

240 

200 

245 

1? 

40 

42 

53 

PEARSON 

10? 

465 

478 

626 

VS 

5? 

258 

242 

302 

N  (Ctl) 

1? 

46 

59 

71 

KOLMOG 

-SHIR  TEST 

20  PTS 

3 C  PTS 

5  0  PTS 

10? 

528 

478 

525 

CONTROL 

5? 

275 

242 

283 

1? 

45 

45 

54 

FEAR  SCN 

10? 

623 

571 

673 

VS 

5? 

339 

310 

371 

N  (C,l) 

1? 

66 

70 

5  5 

FOUTZ 

FN  TEST 

20  PTS 

3C  °TS 

50  FTS 

1C? 

512 

45  8 

5  3  C 

CONTROL 

5? 

248 

207 

263 

1? 

56 

48 

56 

PEARSON 

1C? 

481 

428 

457 

VS 

5? 

241 

196 

234 

N  (0,1) 

1? 

60 

36 

57 

•CONTROL*  IS  IKE  TEST  OF  UNIFORM<0, 1)  VS  UNIFORMS, 1) 

•PEARSON'  IS  THE  PEARSON  DISTRIBUTION 
BETA  1  =  0.0  ,  BETA  2  =  2.30 


MMBEP  OF  HYP  C  THE  S  I  S  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SOURED  TEST 

20  PTS  3 C  PTS  50  PTS 


10? 

435 

443 

541 

CCNTFCL 

5? 

237 

241 

247 

1? 

45 

52 

45 

PEAR  SCN 

10? 

432 

445 

533 

VS 

5? 

243 

247 

241 

N  (C,I) 

1? 

43 

55 

52 

KOLNCG-SNI R 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

1C? 

489 

528 

489 

CCNTFCL 

5 % 

238 

252 

245 

1% 

55 

50 

64 

FEAPSCN 

10 % 

456 

543 

503 

VS 

5% 

242 

274 

262 

N  (C,l) 

1? 

61 

54 

61 

FOUTZ  FN  TEST 

20  PTS  3 C  PTS  50  PTS 


10? 

536 

467 

487 

CCNTFCL 

5% 

275 

22  7 

240 

1? 

61 

34 

43 

FEAPSCN 

10? 

523 

441 

467 

VS 

5% 

265 

217 

228 

N  (C,l) 

1% 

58 

33 

43 

•CCNTFCL’  IS  The  TEST  OF  UMpCRM(OtI)  VS  UNIFORM(Otl) 


'PEARSCN'  IS  THE  PEARSON  DISTRIBUTION 
BETA  1  =  0.0  »  BETA  2  =  2.80 


M.NBEP  CF  HVPC THESIS  REJECTIONS  IN  50C0  REPLICATIONS 
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CHI  SQUARED  TEST 


20  FTS 

3 C  pT S 

50  PTS 

1  c-s 

4C8 

390 

453 

CCNTFCL 

5? 

230 

214 

213 

1 % 

35 

50 

52 

FEAPSCN 

less 

756 

992 

1576 

VS 

5? 

495 

634 

973 

N  <C,1) 

1? 

121 

185 

343 

KOLNCG-SMI  R 

TEST 

20  PTS 

3  C  pTS 

50  FTS 

10? 

492 

515 

520 

CCNTFCL 

5? 

252 

255 

264 

1? 

55 

49 

52 

FEAPSCN 

10? 

942 

1161 

1498 

VS 

5? 

5  46 

693 

944 

N  <C,1) 

1? 

158 

192 

294 

FOUTZ  FN 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

10? 

5  05 

478 

5C4 

CCNTFCL 

5? 

2  23 

261 

2  e  c 

1? 

50 

56 

48 

FEAPSCN 

10? 

8  53 

1400 

2487 

VS 

5? 

469 

884 

1785 

N  (0,1) 

1? 

129 

311 

799 

•CCNTFCL  *  IS  THE  TEST  CF  UNIFCRM(0,1)  VS  U-NI  FORM  (  0 , 1  ) 

•PEARSCN*  IS  THE  PEARSON  DISTRIBUTION 
BETA  1  =  O.Ol,  BETA  2  =  1.75 


NLME  E  F  CF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 

20  PTS  3C  °T S  50  PTS 


10? 

428 

410 

452 

CCNTFGL 

5? 

239 

205 

2  a 

1? 

40 

43 

52 

FEAFSCN 

10? 

614 

744 

1082 

VS 

5? 

375 

43  5 

610 

N  (C,l) 

1? 

93 

120 

153 

KOLMOG-SNIR 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

10? 

543 

470 

5C6 

CCNTFGL 

5? 

250 

229 

274 

1? 

41 

50 

48 

FEAFSCN 

10? 

837 

857 

1106 

VS 

5? 

482 

479 

64  8 

N  <C,1) 

1? 

121 

131 

184 

FOUTZ  FN  TEST 


20  PTS  3C  PTS  50  PTS 


1C? 

5  11 

474 

515 

CCNTFCL 

5? 

240 

214 

259 

1? 

55 

49 

54 

FEAR  SON 

1C? 

558 

726 

1288 

VS 

5? 

305 

405 

774 

N  (C,l) 

1? 

82 

1C1 

233 

•CCNTFCL'  IS  THE  TEST  QP  UN  I PQRM( C , 1 )  VS  UNIFCRM<0,1) 

•PEARSCK'  IS  th£  PEAFSON  DISTRIBUTION 
BETA  1  =  0.01,  BETA  2  =  1.90 


M.H6EF  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATICNS 


86 


CHI  SQUARED  TEST 


20  PTS  3  C  PTS  50  PTS 


CCMRCl 

10? 

403 

442 

510 

5? 

230 

237 

24! 

1? 

48 

41 

47 

FEAR  SCN 

10% 

458 

556 

702 

vs 

5% 

267 

321 

366 

N  (C,I) 

1% 

51 

74 

56 

KOLHOG-SVIR 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

10% 

497 

485 

5C2 

CCNTFCL 

5% 

240 

272 

258 

1% 

55 

50 

65 

PEAR  SCN 

10% 

581 

650 

715 

VS 

5% 

3C9 

2  29 

365 

N  (C,i) 

1% 

74 

73 

54 

FOUTZ  FN  TEST 


20  FT S  30  PTS  50  PTS 


10% 

525 

479 

546 

CCNTFCL 

5% 

241 

241 

285 

1% 

62 

56 

54 

PEARSCN 

10% 

534 

554 

666 

V  S 

5% 

265 

283 

353 

N  (C,l) 

1% 

68 

66 

81 

•CCNTFCL'  IS  THE  TEST  CF  UNIF0RM(0,i)  VS  UNIFORMS, 1) 

•P EARSCN '  IS  THE  PEARSON  DISTRI8UTI0N 
BETA  1  =  0.25,  BETA  2  =  3.20 


M.MEER  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 

20  PTS  3  C  PTS  50  PTS 


10? 

373 

419 

502 

CONTROL 

5? 

22C 

230 

211 

I % 

55 

49 

24 

FEAPSCN 

10? 

739 

1C30 

1672 

VS 

5% 

476 

664 

1099 

N  (C,l) 

\% 

118 

231 

421 

KOLHOG-SNIR 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

10 % 

5  16 

483 

475 

CCNTFCL 

5% 

253 

227 

221 

1 % 

48 

51 

38 

FEAPSCN 

10? 

790 

896 

1125 

VS 

5% 

474 

528 

735 

b  <C,1  > 

1% 

121 

147 

224 

FOUTZ  FN  TEST 

20  PTS  3C  PTS  50  PT$ 


1C? 

499 

481 

477 

CCNTFCL 

5? 

227 

238 

252 

1? 

46 

49 

5C 

FEAPSCN 

10? 

868 

1090 

1516 

VS 

5? 

471 

660 

965 

n  icai 

1? 

121 

185 

344 

‘CCNTFCL'  IS  THE  TEST  OF  UNI  FQRM{  C  ♦  1)  .  VS  UNIFORMS, I) 

' PEARSCN •  IS  THE  PEARSON  DISTRIBUTION 
EET A  1  =  0 .5 0 1  SETA  2  =  3.00 


NLNSEF  OF  HYPOTHESIS  REJECTIONS  IN  50C0  REPLICATIONS 
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CHI  SQUARED 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

1C? 

4  C 1 

428 

522 

CCNTFCL 

5? 

236 

225 

239 

1? 

43 

39 

54 

FEAPSCN 

10? 

2440 

3386 

46  54 

vs 

5? 

1960 

2820 

43  2  C 

N  (C,l) 

1? 

931 

1632 

3267 

KOLMOG-SMIR 

TEST 

20  PTS 

3C  PTS 

50  FT 

1C? 

482 

513 

509 

CCNTFCL 

5? 

239 

254 

245 

1? 

47 

52 

61 

FEAPSCN 

1C? 

11  55 

1459 

20C  1 

VS 

5? 

730 

920 

1428 

N  <  C,1  ) 

1? 

216 

322 

6C9 

FQUTZ  F N  TEST 

20  FT S  30  PTS  50  PTS 


10? 

531 

492 

55C 

CCNTFCL 

5? 

250 

250 

277 

1? 

57 

51 

61 

FEAPSCN 

1C? 

2135 

2837 

392C 

VS 

5? 

1460 

2156 

3309 

N  <  C  ,1  ) 

1? 

583 

1049 

21C4 

'CCNTFCL'  IS  THE  TEST  CF  UNIFORMIO,!)  VS  UNIFORMIO, 1) 

' P  E  ARSCN  '  IS  THE  PEAFSON  DISTRIBUTION 
BETA  1  =  l.OC,  BETA  2  =  3.40 


NLNE  £  F  CF  HVPCTI-ESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SCUiREO  TEST 

2C  PTS  3 C  PTS  50  PTS 


10? 

411 

427 

475 

CCNTFCL 

5? 

226 

238 

227 

1? 

34 

49 

34 

FEARSCN 

1C? 

1446 

2151 

4512 

\  S 

5? 

1070 

2546 

4035 

N  (C,l) 

1? 

4C3 

1368 

2895 

KOL^OG-SMIR 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

10? 

495 

495 

5C4 

CCNTFCL 

5? 

249 

260 

265 

1? 

55 

64 

6  C 

FEARSCN 

1C? 

1032 

1314 

1765 

VS 

5? 

645 

832 

1227 

N  (C,l) 

1? 

211 

265 

453 

FCUTZ  FN 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

10? 

482 

481 

485 

CCNTFCL 

5? 

251 

251 

2C7 

1? 

51 

45 

43 

FEAFSCN 

10? 

1578 

2161 

3115 

vs 

5? 

1009 

1493 

2354 

N  (C,l) 

1? 

348 

633 

1254 

•CCNTFCL'  IS  The  TEST  CP  UNIF0RM(0,1>  VS  UNIFORMS, 1) 

•  PEAR  SCN  '  IS  THE  PEARSON  DISTRIBUTION 
BET*  1  =  1.00,  BETA  2  =  3.60 


NLPBEF  OF  HYPCThESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CHI  SQUARED  TEST 


20  FT  S 

3  C  PTS 

50  PTS 

1C? 

397 

427 

4  6  5 

CCNTFCL 

5? 

2  32 

233 

224 

1? 

41 

60 

56 

FEAPSCN 

1C? 

1086 

2009 

3226 

VS 

5? 

764 

1463 

2542 

N  (C  ,1 J 

1? 

255 

651 

1362 

KCLNCG-SMI R 

TEST 

20  PTS 

3  C  PTS 

50  PTS 

10? 

524 

501 

459 

CCNTFCL 

5? 

267 

250 

252 

1? 

51 

56 

62 

FE/FSCN 

10? 

911 

1150 

1526 

VS 

5? 

5  46 

732 

1022 

N  (Ctl) 

1? 

171 

229 

355 

F0UT2  FN 

TEST 

20  PTS 

3 C  PTS 

50  PTS 

1C? 

521 

460 

525 

CCNTFCL 

5? 

’53 

248 

260 

1? 

50 

63 

46 

FEAP  SCN 

10? 

1261 

1656 

2425 

VS 

5? 

764 

1070 

17C6 

N  (C.l) 

1? 

243 

387 

725 

'CCNTFCL'  IS  TEE  TEST  CF  UMFCRM(0,I)  VS  UNIFCRN<0»1> 

* PEARSCN  *  IS  THE  PEARSON  DISTRIBUTION 
BETA  I  =  l.OC,  BETA  2  =  3. SO 


MReEF  CF  HYPOTHESIS  REJECTIONS  IN  50CO  REPLICATIONS 
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CHI  SQUARED  TFST 

20  PTS  30  PTS  50  PTS 


10* 

371 

431 

464 

CONTROL 

5% 

206 

237 

227 

1% 

27 

48 

45 

FEAPSCN 

10* 

2231 

3C39 

4476 

VS 

5% 

1745 

2282 

3538 

N  i  C  » 1  ) 

1 % 

763 

1242 

2652 

KOLMCG-SMIR 

TEST 

20  PTS 

3 C  DTS 

50  PTS 

1C* 

498 

482 

517 

CCNTFOL 

5* 

250 

261 

261 

1* 

42 

48 

54 

FEAPSCN 

10* 

11  53 

1481 

2082 

VS 

5* 

753 

1005 

1447 

n  (cai 

1? 

221 

353 

621 

FOUTZ  FN  TEST 

20  PTS  30  DTS  50  PTS 


10* 

45C 

498 

486 

CONTROL 

5* 

258 

255 

248 

1* 

44 

45 

61 

FEAPSCN 

1C* 

1872 

2587 

3548 

VS 

5* 

1220 

1920 

2894 

N  <C  ,1  ) 

l* 

469 

886 

1644 

*CCNT  PCI’  IS  THE  TEST  CF  UNIFORM(Otl)  VS  UNIFORMS, 1) 

•FEARSCN*  IS  THE  PEARSON  DISTRIBUTION 
BETA  1  =  2.0C,  BETA  2  =  5.50 


NLMEEP  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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